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Photosynthesis  is  the  major  natural  process  that  provides  free  energy  to  the  metabolism  of 
plants  and  photosynthetic  bacteria.  However,  understanding  of  the  primary  events  in  this 
process  is  still  incomplete.  The  energy  of  incident  photons  is  collected  and  channeled 
into  the  membrane  spanning  reaction  center  where  it  is  used  to  initiate  an  electron 
transfer.  The  first  ion  formed  within  the  bacterial  reaction  center  is  a dimer  of 
bacteriochlorophyll  a molecules,  known  as  the  primary  donor  (D),  which  becomes 
cationic  upon  losing  an  electron.  This  electron  then  begins  its  journey  through  the  A- 
branch  of  the  reaction  center,  which  is  composed  of  a bacteriochlorophyll  molecule,  a 
bacteriopheophytin  molecule,  and  a ubiquinone  molecule.  Upon  leaving  the  primary 
donor,  the  electron  reduces  the  bacteriochlorophyll  a (BChl  a)  molecule  and  then  the 
bacteriopheophytin  a (BPheo  a)  molecule.  The  charge  separation  between  the  primary 
donor  and  the  bacteriopheophytin  has  been  generated  and  studied  by  others.  The  role  of 
BChl  a in  this  electron  transfer  remains  unclear.  In  this  work,  the  radical  cationic  and 
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anionic  states  of  BChl  a and  BPheo  a were  investigated  with  EPR,  ENDOR  and  ESEEM. 
Photosynthetic  plant  pigments,  chlorophyll  a (Chi  a)  and  chlorophyll  b (Chi  b ),  were 
reduced  and  studied  with  high  field  EPR.  This  work  reports  the  first  high  field  EPR  data 
with  fully  resolved  g-tensors  for  the  radical  anions  of  BChl  a,  BPheo  a,  Chi  a,  and  Chi  b. 
Also,  for  the  first  time,  orientationally  selected  ENDOR  was  used  to  measure  the  proton 
hyperfine  interactions  of  BChl  a and  BPheo  a radical  anions.  These  anions  show  similar 
trends  in  their  anisotropic  hyperfine  interactions,  but  the  differences  in  the  Zeeman  and 
hyperfine  interactions  are  discemable  at  high  fields.  The  orientations  and  values  of  the 
Zeeman  and  hyperfine  tensors  have  been  approximated  using  spectral  fitting  procedures 
and  density  functional  theory  predictions.  This  work  marks  the  first  steps  toward  full 
characterization  of  the  electronic  nature  of  the  primary  acceptors  in  (bacterial  and  plant) 
photosynthesis. 
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CHAPTER  1 
PHOTOSYNTHESIS 


Photosynthesis  is  the  process  used  by  plants,  algae,  and  photosynthetic  bacteria  to 
convert  the  energy  of  light  into  chemical  energy.  This  energy  is  used  in  different  ways 
depending  on  the  specific  organism.  Green  plants  and  algae  (eukaryotes)  use 
photosynthesis  to  consume  carbon  dioxide  and  produce  oxygen  and  carbohydrates,  while 
certain  bacteria  (prokaryotes)  consume  organic  compounds  and  can  adapt  to  their 
environment  to  survive  in  a variety  of  conditions.  Basically,  photosynthesis  provides  the 
energy  and  reduced  carbon  necessary  for  life  on  this  planet,  as  well  as  the  oxygen 
necessary  for  the  survival  of  aerobic  organisms.  The  last  200  years  have  led  to  great 
achievements  in  understanding  the  basic  processes  of  photosynthesis,  yet  fundamental 
uncertainties  remain  [1-6].  After  a brief  historical  perspective,  the  structure  of  the 
photosynthetic  machinery  will  be  discussed  and  then  a closer  look  will  be  taken  at  the 
chemistry  of  this  process. 

The  first  step  toward  understanding  the  process  of  photosynthesis  came  in  the 
discovery  of  oxygen  production.  In  1780,  Joseph  Priestly  kept  a mouse  alive  through  the 
oxygen  production  of  a sprig  of  mint  [7].  Priestly  also  showed  that  a burning  candle 
would  go  out  when  under  a closed  vessel,  but  the  presence  of  a sprig  of  mint  would  allow 
the  candle  to  bum  again  after  several  days.  He  did  not  know  about  molecular  oxygen,  but 
he  did  ultimately  claim  that  plants  could  “restore  air  which  had  been  injured  by  the 
burning  of  candles”  [8,  p.  432],  By  the  beginning  of  the  19th  century,  understanding  of 
photosynthesis  was  sufficient  to  develop  the  basic  equation  we  know  today. 
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C02  + 2H2X  + hv  = (CH20)  + X2  + H20 

The  energy  needed  to  drive  this  reaction  is  supplied  by  sunlight.  The  X in  this  reaction 
denotes  oxygen  for  plant  photosynthesis  and  sulfur  for  certain  photosynthetic  bacteria, 
such  as  Rb.  sphaeroides  and  Rps.  viridis  [7].  This  rather  simple  equation  represents  an 
amazing  orchestration  of  macromolecules  and  proteins,  which  carry  out  a very  efficient 
conversion  of  solar  energy  to  chemical  energy.  In  1845,  Julius  Robert  von  Mayer 
proposed  that  photosynthetic  organisms  convert  light  energy  into  chemical  free  energy 

[8] .  The  magnitude  of  Mayer’s  statement  is  still  being  realized  today,  since  the  fine 
details  of  this  energy  conversion  process  are  still  not  fully  understood.  In  the  1930’s, 
Gaffron  and  Wohl  postulated  that  incident  light  is  absorbed  by  antenna  chlorophyll 
molecules  to  produce  an  exciton,  which  migrates  to  a more  stable  photochemical  trapping 
chlorophyll  molecule,  which  is  part  of  the  integral  photosynthetic  unit:  the  reaction  center 

[9] ,  It  was  not  until  the  1980’s  that  the  structures  of  bacterial  reaction  centers  were 
solved  [6,10-13],  and  it  was  2000  before  plant  photosystem  crystal  structures  were  solved 
[14,15].  Solving  the  structure  of  the  reaction  center  provided  a spatial  reference  of  the 
components  involved  in  the  energy  conversion,  which  allowed  research  to  focus  on  the 
energetics  and  electronic  structure. 

In  1956,  Commoner  et  al.  measured  a radical  species  with  EPR  generated  from 
photo-oxidized  chloroplasts  [16].  The  nature  of  this  ionic  species  would  be  debated  for 
years,  until  it  was  identified  as  a chlorophyll  cation  radical  [17,18-21,22-25],  In  the  early 
1970’s,  this  species  was  further  identified  as  a dimer  of  bacteriochlorophyll  molecules  in 
bacterial  reaction  centers  through  EPR  and  ENDOR  studies  [26-29,30-32],  as  well  as 
through  X-ray  crystallography  [6,33],  Further  magnetic  resonance  studies  have  examined 
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the  electronic  structure  of  this  dimer  to  estimate  the  properties  and  location  of  the 
electron  density  [34,35]. 

Energy  Conversion 

Photosynthesis  can  be  divided  into  three  separate  but  complementary  processes: 
collection  of  incident  light  and  transfer  of  this  light  energy  to  special  trans-membrane 
complexes  known  as  reaction  centers,  charge  separation  within  the  reaction  center 
through  electron  transfer,  and  protonation  of  electron  acceptors  and  establishment  of  a 
proton  gradient  [36].  The  work  presented  here  will  focus  on  the  charge  separation  and 
electron  transfer  pathway.  To  appreciate  the  pathway  of  the  electron  transfer  reactions, 
the  structure  of  the  photochemical  machinery  must  be  understood. 

Solar  energy  enters  the  reaction  center  through  a pair  of  chlorophyll  molecules 
known  as  the  special  pair  or  the  primary  donor.  The  primary  donor  (D)  uses  this  excess 
energy  to  initiate  an  electron  transfer  within  a chain  of  cofactor  molecules,  which 
comprise  the  heart  of  the  reaction  center  (RC).  This  RC  complex  is  composed  of  cofactor 
molecules  and  proteins,  which  perform  an  efficient  electron  transfer  through  the  stromal 
(for  green  plants)  or  cytoplasmic  (for  bacteria)  membrane.  Electron  transfer  reactions 
occur  between  the  primary  donor  and  intermediate  acceptor  molecules,  creating  charge 
separated  pairs  between  the  primary  donor  and  each  consecutive  intermediate  acceptor. 
The  ultimate  charge  separation  is  equalized  by  a proton  gradient  to  drive  the  production 
of  ATP.  This  series  of  oxidative  and  reductive  steps  is  heavily  driven  by  the  redox 
potential  and  molecular  orbital  overlap  of  each  consecutive  cofactor.  Although  the 
mechanism  of  the  photosynthetic  process  differs  between  green  plants  and  bacteria,  the 
efficient  conversion  of  light  into  chemical  energy  is  similar.  Photosynthetic  bacteria 
contain  only  one  type  of  reaction  center,  while  green  plants  have  two:  photosystem  I and 
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photosystem  II  [37],  These  photosystems  work  in  series  with  photosystem  I being  the 
primary  energy  source  for  the  plant  while  photosystem  II  is  mainly  responsible  for 
producing  oxygen  and  splitting  water  [37],  Since  photosystem  I has  a strong  optical 
absorbance  at  X ~ 700  nm,  and  photosystem  II  requires  light  of  X < 680  nm,  the  primary 
donor  molecule  that  initiates  the  electron  transfer  within  each  photosystem  is  termed  P700 
and  P6so,  respectively.  Bacterial  primary  donors  are  labeled  similarly  as  P865  for 
Rhodobacter  sphaeroides  and  P960  for  Rhodopseudomonas  viridis.  Bacterial  reaction 
centers  build  the  primary  donor  from  bacteriochlorophyll  a and  b molecules  for  Rb. 
sphaeroides  and  Rps.  viridis,  respectively.  These  chlorophyll  molecules  are 
photochemical  traps  for  the  excess  energy  of  the  light  harvesting  system.  The  energy 
conversion  in  the  reaction  center  entails  the  absorption  of  one  photon,  which  results  in  a 
near  unity  quantum  yield  [6],  This  efficiency  also  applies  to  photosynthetic  bacteria, 
which  seem  to  have  certain  structural  and  functional  similarities  to  photosystem  II 
[38,39].  The  crystal  structure  of  plant  photosystems  has  just  recently  been  solved 
[14,15,37],  which  makes  it  advantageous  to  turn  to  certain  non-sulfur  bacterial  reaction 
centers  whose  reaction  center  structures  are  known  to  high  resolution,  such  as  Rb. 
sphaeroides  [13]  and  Rps.  viridis  [6,10-12,40], 

The  Bacterial  Reaction  Center 

The  bacterial  reaction  center  is  composed  of  a protein  scaffolding  and  cofactor 
molecules  which  perform  the  chemistry  of  the  energy  conversion.  The  cofactor 
components  of  the  reaction  center  are:  a bacteriochlorophyll  dimer,  two 
bacteriochlorophyll  monomers,  two  bacteriopheophytin  molecules,  two  quinone 


5 


Figure  1-1  Rb.  sphaeroides  reaction  center  crystal  structure  [41]  consisting  of  a primary 
donor  (D),  a carotenoid  (Car),  two  accessory  bacteriochlorophyll  molecules 
(B),  two  bacteriopheophytin  molecules  (®),  two  ubiquinone  molecules  (Q), 
and  a non-heme  iron  (Fe). 

molecules,  and  a non-heme  iron  atom.  Figure  1-1  is  a crystal  structure  image  of  the 
bacterial  reaction  center  from  Rb.  sphaeroides  [13].  Branches  A and  B are  chemically 
identical,  with  the  exception  of  a carotenoid  molecule  on  the  B branch  which  quenches 
the  triplet  state  on  the  primary  donor  (3D*)  [42].  These  branches  maintain  an 
approximate  C2  symmetry  about  an  axis  running  between  the  branches,  through  the 
planes  of  the  primary  donor  and  the  Fe  atom  [10,43,44],  Despite  structural  symmetry  in 
the  pigments  of  the  reaction  center,  the  electron  transfer  occurs  asymmetrically, 
following  only  the  A branch  [45,46].  The  origins  of  this  behavior  have  been  postulated 
to  result  from  the  alteration  of  the  MO  on  bacteriopheophytin  (®  A)  by  the  protein 
environment  [47,48].  The  LUMO  of  <t>A  must  overlap  with  the  HOMO  of  BA\  in  order 
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for  an  electron  transfer  to  occur  between  these  two  molecules  [49].  Electrostatic 
influences,  such  as  polar  protein  side -chains,  are  suspected  to  play  a role  in  the 
unidirectionality  of  the  electron  transfer.  Depending  on  the  rotation  angle  of  the  acetyl 
group  on  ring  I,  hydrogen  bonding  interactions  can  affect  the  LUMO,  which  affects  the 
electron  transfer  [49-51],  Upon  receiving  the  solar  energy  from  the  antenna  chlorophylls, 
the  primary  donor  (D)  promotes  an  electron  within  the  singlet  manifold  to  the  first 
excited  state.  Upon  leaving  the  primary  donor,  the  electron  density  generates  an  anionic 
state  on  a bacteriochlorophyll  molecule  (BA),  a bacteriopheophytin  molecule  (Oa),  and  a 
quinone  molecule  (QA). 

This  electron  transfer  creates  a charge  separated  state  between  D+  and  possibly 
Ba'  [4,52-59]  which  very  quickly  results  in  the  reduction  of  Oa  to  create  D+ — <1>A‘  [53- 
55,60].  The  exact  role  of  BA  in  this  process  is  still  under  discussion.  The  initial  charge 
separation  between  D and  BA  occurs  in  about  3ps  [61-64],  This  charge  separated  state, 
D+ — Ba',  state  decays  in  less  than  1 ps  to  generate,  D+ — <J>A‘  [53,54,65,66].  The 
electron  transfer  then  proceeds  to  QA  in  about  200  ps  [6,41].  The  rate  of  these  reactions 
indicates  the  efficiency  of  this  process,  but  introduces  an  element  of  difficulty  in  the 
study  of  this  process.  Figure  1-2  illustrates  the  relative  energies  of  the  various  ionic 
species  involved  in  the  electron  transfer  process.  The  approach  of  this  work  will  be  to 
isolate  various  ionic  species  to  study  the  electronic  properties  of  the  individual  radicals, 
which  dictates  the  chemistry  of  this  electron  transfer  process. 
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Figure  1-2  Relative  energy  levels  for  the  radical  pairs  in  the  primary  photochemistry  in 
photosynthetic  bacteria.  Relative  energies  are  not  to  scale. 

Ionic  Cofactors 

The  reaction  center  cofactors  have  slight  molecular  variations,  which  are  specific  to 
the  organism  in  which  they  are  found.  Rb.  sphaeroides  is  composed  of  a-type 
bacteriochlorin  pigments,  Rps.  viridis  is  composed  of  fr-type  bacteriochlorin  pigments, 
and  eukaryotic  reaction  centers  are  composed  of  a-type  and  b-type  chlorin  pigments 
(shown  in  Figure  1-3).  The  main  chemical  difference  between  prokaryotic 
(bacteriochlorin)  and  eukaryotic  (chlorin)  pigments  lies  in  the  extent  of  saturation  on  the 
tetrapyrrole  [67].  Bacteriochlorin  molecules  have  two  saturated  pyrrole  rings  (rings  II 
and  IV),  while  chlorin  molecules  have  one  reduced  pyrrole  ring  (ring  IV).  This 
difference  reduces  the  extent  of  conjugation  on  the  plant  pigments,  which  increases  the 
oxidation  and  reduction  potential  for  the  plant  species,  relative  to  bacterial  pigments. 
Within  the  reaction  centers,  the  protein  environment  influences  the  redox  behavior  as 
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well,  but  the  nature  of  the  individual  species  will  be  the  focus  here.  According  to  Marcus 
theory,  the  rate  and  efficiency  of  the  radical  pairing  depend  on  the  extent  of  molecular 
orbital  (MO)  overlap,  the  difference  in  redox  energy,  and  the  reorganization  energy  [68]. 
The  cationic  state  serves  as  a direct  probe  of  the  highest  occupied  molecular  orbital 
(HOMO),  while  the  anionic  state  probes  the  lowest  unoccupied  molecular  orbital 
(LUMO).  Hence,  the  intermediate  acceptor  radical  cation  and  anion  species  have  been 
isolated  from  the  reaction  center  and  studied  with  electron  magnetic  resonance  (EMR). 

In  this  way,  the  frontier  orbitals  may  be  studied,  which  are  responsible  for  the  electron 
transfer  pathway. 

The  recent  crystallization  of  photosystems  I and  II  have  identified  the  cofactors 
responsible  for  the  primary  charge  separation  in  higher  plants  [14].  The  apparent 
versatility  of  chlorophyll  a in  these  photosystems  prompted  the  study  of  this  species.  Chi 
a seems  to  be  a plant  analogue  to  BChl  a in  bacterial  reaction  centers.  The  absence  of 
chlorophyll  b in  the  photosystems  generated  curiosity,  and  applied  magnetic  resonance  to 
probe  any  differences  in  electronic  structure  between  these  anions. 
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Bacteriochlorophyll  a 


Bacteriopheophytin  a 


Chlorophyll  a Chlorophyll  b 

Figure  1-3  Molecular  structures  of  prokaryotic  (bacterial)  and  eukaryotic  (plant) 
pigments. 


CHAPTER  2 

MAGNETIC  RESONANCE  THEORY 


Magnetic  resonance  spectroscopy  is  a technique  that  uses  an  external  magnetic 
field  to  lift  the  degeneracy  of  spatial  quantization  among  the  dipoles  in  the  sample. 
Nuclear  Magnetic  Resonance  (NMR)  is  a popular  technique,  which  probes  the  nuclear 
magnetic  moments  of  a system,  while  Electron  Magnetic  Resonance  (EMR),  also  known 
as  Electron  Spin  Resonance  (ESR)  or  Electron  Paramagnetic  Resonance  (EPR),  probes 
the  electronic  dipole  moments  of  a system.  Since  the  EMR  performed  in  this  work, 
studies  paramagnetic  radical  ions,  the  expression  EPR  will  be  applied  from  now  on. 
Certain  nuclei  are  magnetic,  and  thus,  can  be  investigated  using  NMR.  Magnetic  nuclei 
are  commonly  found  in  organic  molecules,  such  as  hydrogen  and  nitrogen,  which  gives 
NMR  a certain  utility  among  synthetic  organic  chemists.  All  electrons  on  the  other  hand, 
possess  a magnetic  moment.  The  interaction  of  these  nuclear  and  electronic  moments 
contributes  to  the  utility  of  EPR.  Electrons  occupy  orbitals,  which  are  filled  according  to 
Hund’s  rules.  These  rules  require  electron  pairing  within  orbital  sublevels,  which  leads 
to  a cancellation  of  the  total  spin.  If  at  least  one  electron  remains  unpaired,  the  molecule 
is  paramagnetic  and  may  be  investigated  with  EPR  to  study  molecular  structure  and 
geometry,  inter-molecular  interactions,  and  molecular  orbital  energy  level  behavior.  The 
main  focus  here  will  be  the  application  of  EPR  to  photosynthetic  pigment  ions  to 
determine  structural  information  and  predict  chemistry  of  the  various  ionic  species 
involved  in  the  primary  charge  separation  of  bacterial  and  plant  reaction  centers.  A 
variety  of  EPR  experiments  were  used  to  investigate  these  systems,  such  as,  Electron 
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Spin  Echo  Envelope  Modulation  (ESEEM),  high  field  cw-EPR,  and  Electron  Nuclear 
DOuble  Resonance  (ENDOR).  An  introduction  to  EPR  will  be  limited  but  a myriad  of 
texts  exists  for  further  background  information  [69-71]. 

Magnetic  Moments  and  Angular  Momenta 
Electrons  are  magnetic  particles  that  will  interact  with  an  externally  applied 
magnetic  field.  Classical  physics  predicts  the  moment  will  precess  about  an  externally 
applied  magnetic  field  with  a certain  angular  velocity.  This  classical  picture  would  allow 
the  moment  to  maintain  any  arbitrary  angle  with  the  external  field  during  its  precession. 

In  the  quantum  picture,  however,  the  angle  between  the  moment  and  the  external  field, 

B0,  is  quantized,  which  only  allows  space-quantized  values  of  angular  momentum.  Since 
the  electron’s  magnetic  moment  and  angular  momentum  are  directly  related  to  the 
energy,  this  quantization  of  momentum  leads  to  a quantization  of  energy.  Figure  2-1 
illustrates  the  effect  of  a strong  magnetic  field  on  the  electronic  dipoles  within  a sample. 
When  B0  is  turned  on,  the  electrons  align  with  this  field  and  precess  about  the  axis  of  B0 


B0 


Figure  2-1  The  effect  of  a strong  external  magnetic  field,  B0,  on  dipole  moments  in  a 
paramagnetic  sample.  The  vectors  in  the  left  frame  represent  individual 
dipoles,  while  the  larger  vectors  in  the  right  frame  are  net  dipoles,  resulting 
from  the  electronic  dipole  precession  about  B0. 
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with  spatially  quantized  orientations.  The  vectors  shown  on  the  right  of  figure  2-1  are 
actually  resultants  of  each  electronic  dipole  moment’s  precession  about  B0.  The  arrows 
parallel  to  B0  represent  dipoles  in  the  lower  spin  state,  while  the  antiparallel  arrows 
represent  dipole  moments  in  the  upper  energy  state.  The  frequencies  of  precession  will 
be  intimately  related  to  the  nature  of  the  electron  density  on  the  paramagnetic  molecules 
in  the  sample.  This  relationship  is  described  through  the  electronic  Zeeman  equation. 

HZeemtm=\iB-g-B0-S  (2.1) 

/iv  = pB-g-B0  (2.2) 

Equation  2.1  is  the  Zeeman  component  of  the  Hamiltonian  equation,  where  equation  2.2 
describes  the  Zeeman  resonance  condition,  or  difference  between  two  electronic  Zeeman 
levels.  The  left  side  of  this  equation  is  the  microwave  energy  and  the  right  side  is  the 
Zeeman  energy  splitting  that  is  proportional  to  the  external  magnetic  field.  In  this 
equation  h represents  Planck’s  constant  (6.626  x 10  34  Js),  v represents  the  microwave 
frequency  (in  Hz),  g represents  the  g-tensor,  pB  is  the  Bohr  magneton  (9.27402  x 10‘24 
JT1),  S is  the  spin  angular  momentum,  and  B0  is  the  static  external  magnetic  field.  When 
the  Zeeman  condition  is  satisfied,  an  EPR  transition  occurs,  as  shown  in  figure  2-2.  The 
selection  rule  limits  transitions  to  AMs  = ±1,  which  explains  why  S does  not  appear  in 
equation  2.2. 

Electrons  possess  two  forms  of  angular  momentum:  orbital  and  spin.  Thus,  the 
total  angular  momentum  is  given  as  the  sum  of  the  orbital  and  spin  components.  Orbital 
angular  momentum  (L)  can  be  accounted  for  quantum  mechanically  as  well  as 
classically.  In  the  classical  picture,  this  quantity  results  from  the  electron’s  motion  within 
the  orbitals  of  the  molecule.  This  motion  generates  a magnetic  field,  which  alters  the 
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Figure  2-2  (a)  Net  spin  vector  (which  has  a resulting  dipole  in  the  opposite  direction)  and 
(b)  energy  level  scheme  for  S = Vi  system,  having  two  allowed  energy  states. 

magnetic  field  felt  by  neighboring  molecules  in  the  sample.  Spin  angular  momentum  (S), 

on  the  other  hand,  is  purely  a relativistic  quantum  property  that  has  been  detected 

experimentally.  The  coupling  of  the  spin  (S)  and  orbital  (L)  angular  momenta  is 

measured  by  the  value  of  a sample’s  g-tensor.  All  the  species  to  be  dealt  with  here  have 

one  unpaired  electron,  S=  V2,  that  is  delocalized  over  the  conjugated  7t-system  on  the 

molecule,  which  has  relatively  weak  spin-orbit  interactions  leading  to  g = 2.  Spin  angular 

momentum  can  also  interact  with  neighboring  dipoles— nuclear  ( I ) or  electronic  (S)~ 

which  changes  the  local  field  felt  by  the  electron  spin.  This  dipolar  interaction  is  known 

as  hyperfine  (I-S)  or  fine  structure  (S-S)  interactions  and  will  be  discussed  after  EPR. 

cw-EPR 

An  EPR  experiment  may  be  performed  in  one  of  two  ways,  continuous  wave  (cw) 
or  pulsed.  In  the  continuous  wave  experiment  the  sample  is  continuously  irradiated  with 
microwaves  of  a fixed  frequency.  The  magnetic  component  of  this  radiation,  Bi, 
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interacts  with  the  electronic  magnetic  moments  of  the  sample  while  the  static  field,  Bc,  is 
swept  through  the  resonance  condition  described  by  equation  2.2.  Figure  2-2  illustrates 
the  behavior  of  the  net  spin  magnetization  vector  in  the  presence  of  the  magnetic 
components  of  the  external  static  field,  B0,  and  the  microwave  radiation,  Bj.  When  the 
frequency  of  the  incident  microwave  radiation  and  the  static  field  satisfy  the  resonance 
condition,  radiation  is  absorbed  by  the  sample.  This  radiation  is  absorbed  by  the  spins 
that  are  “flipped”  from  the  lower  energy  state  to  a higher  energy  state.  That  is,  when  the 
resonance  condition  is  satisfied  a portion  of  the  spin  population  is  promoted  from  a lower 
energy  state  to  a higher  energy  state,  which  reduces  the  population  of  the  lower  level  and 
increases  the  population  of  the  upper  level.  The  relative  populations  of  these  energy 
states  determine  the  intensity  of  the  microwave  absorption.  If  the  upper  state’s 
population  approaches  the  population  of  the  lower  state,  saturation  effects  distort  the  EPR 
signal  and  can  reduce  its  intensity.  The  number  and  splitting  of  these  energy  levels 
depend  on  the  number  and  nature  (orbital  location)  of  the  unpaired  electron(s),  and  on  the 
magnetic  nuclei  interacting  with  the  electron  density.  These  electronic  energy  states  may 
interact  with  neighboring  magnetic  nuclei  to  produce  hyperfine  splitting  of  the  electronic 
Zeeman  states. 

g-Tensor 

The  g-tensor  is  the  proportionality  constant  that  represents  the  Zeeman  levels’ 
deviation  from  a true  spin  system.  This  quantity  can  be  single  valued  (isotropic),  or  can 
be  anisotropic  due  to  anisotropic  spin-orbit  interactions  and  mixing  of  the  ground  state 
with  excited  states.  The  g-tensor  can  be  a very  useful  parameter  for  identifying 
molecules,  determining  geometrical  molecular  structure,  or  assessing  interactions 
between  a radical  and  its  environment  [33,72-75].  The  interaction  of  ground  and  excited 
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electronic  states  also  influence  this  parameter.  The  g-value  for  a “free”  electron,  which  is 
free  of  orbital  interaction,  is  known  very  accurately  to  be:  2.002319304386(20)  [76]. 
Organic  free  radicals  maintain  isotropic  g-values  close  to  the  free  electron  value. 
However,  interaction  with  excited  states  and  spin-orbit  coupling  introduce  slight 


relatively  small.  Equation  2.3  shows  the  relationship  between  the  g-value  and  the  free 
electron: 


This  equation  results  from  Rayleigh-Schroedinger  perturbation  treatment  of  the  spin 
Hamiltonian,  which  represents  the  individual  elements  of  the  g-matrix,  and  their 
dependence  on  the  orbital  functions  of  the  excited  states,  n.  The  resulting  g-value,  gm„  is 
equal  to  the  free  electron  g-value,  ge,  unless  orbital  angular  momentum  is  obtained  from  n 
excited  states,  which  is  represented  by  the  summation  of  excitations  to  higher  states.  The 
denominator  of  the  mixing  term  is  the  energy  difference  between  the  excited  state,  E„  and 
the  ground  state,  E0.  The  spin-orbit  coupling  constant,  X,  is  unique  to  individual  nuclei. 

Initially  proposed  by  Stone  [77]  and  modified  by  Pryce  [78],  it  was  determined 
that  the  component  of  the  g-tensor  normal  to  the  plane  of  conjugation  experiences 
minimal  spin-orbit  coupling,  and  thus,  falls  closest  to  the  free  electron  value.  Anisotropic 
Zeeman  interaction  gives  rise  to  the  characteristic  powder  pattern  shown  in  figure  2-3. 
The  theory  of  this  interaction  has  been  further  developed  for  conjugated,  doublet  71- 
radicals  to  characterize  molecular  structure  [41,75,79-81]  and  intermolecular  interactions 


anisotropy  in  the  g-tensor,  but  this  anisotropy  is  usually  on  the  order  of  10'3,  which  is 


(2.3) 


[74,81-83].  Advances  in  high  magnetic  field  technology  have  turned  the  focus  of  much 
research  to  better  understanding  this  parameter.  In  doublet  radicals  with  few  heavy  atoms 
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and  a largely  (spatially)  delocalized  electronic  wavefunction,  the  g-tensor  is  usually  too 
isotropic  to  detect  with  conventional,  X-band  (9.5  GHz)  experiments,  showing  anisotropy 
on  the  order  of  10"4.  Often,  frequencies  of  at  least  200  GHz  are  required  to  obtain  full 
resolution  in  g.  The  principal  components  of  the  g-tensor  can  be  measured  if  the 
anisotropy  in  the  Zeeman  interaction  exceeds  the  unresolved  hyperfine  (nuclear) 
interaction.  The  g-tensor  can  be  theoretically  predicted  as  a function  of  spin-orbit 
interactions  and  electronic  Zeeman  interactions,  with  relativistic  mass  and  gauge 
corrections  [75,79].  The  interaction  of  the  electronic  spin  density  with  the  magnetic  field 
generated  by  the  electron’s  orbital  motion  gives  rise  to  spin-orbit  interactions.  This 
parameter  will  be  discussed  further  in  the  following  section. 


Hyperfine  interaction  is  expressed  in  the  Hamiltonian  as  the  sum  of  the  products  of 
the  spin  moment  vector,  S,  the  hyperfine  coupling  constant,  A„  and  the  nuclear  moment 
vector,  I„  as  shown  in  the  following  equation. 


This  summation  is  evaluated  over  all  nuclei,  i,  and  represents  the  hyperfine  interaction 
between  an  unpaired  electron  and  magnetic  nuclei.  This  interaction  is  independent  of 
field;  however,  the  nuclear  Zeeman  interaction  is  field  dependent.  In  solution,  the 
hyperfine  interaction  becomes  isotropic  due  to  the  averaging  effect  of  molecular  motion. 
In  the  disordered  solid  phase,  such  as  frozen  solution,  the  anisotropic  hyperfine 
interaction  on  large  molecules  with  many  weak  nuclear  interactions,  such  as  chlorophylls, 
produce  such  weak  coupling  that  the  EPR  spectral  line  becomes  broadened  by  this  effect. 
Chlorophylls  possess  approximately  6 alpha,  15  beta,  and  8 gamma  protons,  as  well  as  4 


Hyperfine  Interaction 


(2.4) 
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nitrogen  atoms  that  can  contribute  to  the  hyperfine  interaction  (see  figure  1-3).  Since  the 
electron  density  is  delocalized  over  the  n system  of  about  20  carbon  atoms,  the  spin 
density  on  each  atom  is  relatively  small,  leading  to  relatively  small  hyperfine  couplings 
from  the  aforementioned  nuclei.  In  addition,  the  coupled  nuclei  are  not  symmetry 
related,  and  thus  maintain  slightly  different  hyperfine  coupling  constants.  So,  the  number 
and  diversity  of  hyperfine  interactions,  in  addition  to  the  small  g-anisotropy,  lead  to  an 
unresolved  EPR  line  at  X-band  with  a linewidth  -12-18  G.  This  lack  of  resolution 
prevents  accurate  g-tensor  measurements  and  must  be  overcome  with  higher  fields.  Even 
though  Hhfi  is  independent  of  the  external  field,  B0,  another  term  in  the  Hamiltonian  must 
be  considered.  Consider  the  entire  Hamiltonian  for  this  system: 


^ ^ Zeeman  ^ NuclearZeeman  ^ hfi 

H = m,  g.  B„  S-  -gN.  I,+  A,- 1, 


(2.5) 


If  the  g-tensor  is  not  fully  resolved,  the  inhomogeneous  EPR  linewidth  (ABi/2)  must  be 
dominated  by  field  independent  factors,  such  as  Hhfi.  By  increasing  B0,  the  Hzeeman  term 
and  the  HNuciearzeeman  term  will  eventually  overcome  the  Hhfi  term,  allowing  for  full 
resolution  in  g.  The  magnitude  of  the  field  required  to  fully  resolve  g depends  on  the 
anisotropy  of  the  g-tensor  as  well  as  the  magnitude  of  the  hyperfine  interaction.  For  a 
given  hyperfine  interaction,  the  field  required  to  resolve  the  g-tensor  components 
decreases  as  the  g-anisotropy  increases.  Similarly,  if  the  g-anisotropy  is  very  small,  the 
external  field  required  for  full  resolution  will  be  large.  The  experimental  aspects  of  high 
field  experiments  have  been  thoroughly  reviewed  by  others  [5]  and  will  only  be 
highlighted  here. 

The  spectral  separation  (AB)  in  g-values,  gi  and  g2,  is  given  by  equation  2.6. 
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M-b  v&i  82  j 


(2.6) 


The  ability  to  resolve  the  g-tensor  components  is  enhanced  by  a factor  of  40  when 
increasing  the  microwave  frequency  from  9.5  GHz  to  330  GHz.  G-tensor  resolution  is 
determined  by  the  magnitude  of  the  anisotropic  hyperfine  interactions.  The  magnetic 
nuclei,  consisting  mainly  of  protons  and  nitrogens  in  photosynthetic  pigments,  can 
interact  with  the  magnetic  moment  of  the  electron  density.  This  interaction  splits  the 
electronic  Zeeman  energy  levels  into  hyperfine  levels.  The  transitions  between  these 
levels  can  be  partially  resolved  in  solution  state  EPR  since  the  molecular  tumbling 
averages  out  the  anisotropic  interaction.  If  the  linewidth  of  the  anisotropic  hyperfine 
interaction  exceeds  the  splitting  between  the  g-tensor  components,  then  the  g-tensor  is 
not  resolved  and  higher  fields  must  be  employed. 


Equation  2.7  expresses  the  “high  field”  condition,  which  occurs  when  the  anisotropic 
Zeeman  interaction  exceeds  the  linewidth  of  the  unresolved  hyperfine  interactions.  The 
electron  density  occupies  the  molecular  orbitals  that  participate  in  the  71-conjugation. 
This  MO  has  a similar  orientation  on  each  paramagnetic  molecule  in  the  sample.  In 
randomly  oriented  samples,  such  as  a frozen  solution,  the  molecular  axes  are  distributed 
over  all  possible  orientations  with  respect  to  the  external  field  axis.  Equation  2.8 
illustrates  the  spatial  relationship  between  the  effective  g-value,  g(0,(|> ) and  the  external 
magnetic  field: 


B 


8 iso 


(2.7) 


(2.8) 
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This  angular  distribution  creates  a characteristic  EPR  spectrum  known  as  a powder 
pattern  [63],  Figure  2-3  contains  two  possible  shapes  for  the  g-tensor  to  adopt.  A 
rhombic  geometry  has  three  unequal  diagonal  elements  in  the  tensor,  gxx  > gyy  > gzz,  while 
an  axial  geometry  resembles  a cylindrical  shape  with  two  equal  diagonal  elements,  gxx  = 
gyy>g  ZZ*  The  edges  of  these  absorption  spectra  result  from  distinct  sets  of  molecules 
with  a single  orientation  relative  to  the  external  field.  That  is,  only  molecules  with  the  x- 
axis  of  their  g-tensor  parallel  to  B0  absorb  at  the  gxx  point  of  the  spectrum.  This  condition 
applies  to  the  high  field  side  of  the  spectrum,  where  molecules  with  gzz  parallel  to  B0  are 
in  resonance.  Full  resolution  of  the  g-tensor  is  a goal  of  this  work  and  has  important 
implications  that  will  be  addressed  when  high  field  ENDOR  is  discussed. 


Figure  2-3  EPR  powder  patterns  showing  (a)  rhombic  (gx  > gy  > gz)  and  (b)  axial 
symmetry  (g1>g|| ). 
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cw-ENDOR 

Electron  Nuclear  DOuble  Resonance  (ENDOR)  was  experimentally  discovered  by 
George  Feher  in  1959  [84].  This  experiment  monitors  the  intensity  of  an  EPR  transition 
while  stimulating  NMR  transitions,  so  it  is  basically:  EPR  detected  NMR.  If  hyperfine 
coupling  exists,  the  EPR  signal  intensity  will  be  affected  by  the  NMR  transitions.  This 
experiment  saturates  the  EPR  transition,  then  sweeps  radiofrequency  radiation  to  induce 
NMR  transitions.  If  the  nuclei  experiencing  NMR  are  coupled  to  the  electron  density,  the 
intensity  of  the  EPR  absorption  will  change  as  a function  of  RF  frequency.  Figure  2-4 
depicts  an  S='/2  system  coupled  to  a proton,  I=Vi,  demonstrating  the  theory  of  cw- 
ENDOR.  The  Boltzman  population  is  disturbed  by  microwave  radiation  intense  enough 
to  saturate  (equalize)  the  populations  in  the  upper  and  lower  states.  Since  the  intensity  of 
the  EPR  transition  is  proportional  to  the  difference  in  populations 
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Figure  2-4  cw-ENDOR  population  effects  on  a system  with  S=‘/2,  1=1/2  at  thermal 
equilibrium,  with  relative  spin  population  depicted  as  dots. 
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between  the  states,  any  change  in  the  populations  between  the  energy  states  involved  in 
the  transition  will  affect  the  intensity  of  the  EPR  signal.  Hence,  when  NMR  is  achieved, 
spins  are  taken  out  of  the  upper  level  (or  pumped  into  the  lower  level)  to  enhance  the 
population  difference  between  EPR  levels. 

ENDOR  experiments  can  shed  light  on  the  identity  of  coupled  nuclei,  number  of 
coupled  nuclei,  and  amount  of  spin  density  present  at  each  coupled  nucleus.  In  a powder, 
or  frozen  solution,  the  system  is  randomly  oriented  and  the  ENDOR  resonance  lines  are 
broadened  by  the  different  orientations  between  the  hyperfine  tensors,  A-tensors,  of  the 
individual  nuclei  and  the  external  magnetic  field.  Increasing  the  resolution  of  the  g- 
tensor  can  increase  the  resolution  in  the  ENDOR  spectra  with  orientational  selection. 

This  technique  performs  ENDOR  experiments  at  the  turning  points  in  the  powder  pattern. 
Each  turning  point,  corresponds  to  the  radicals  in  the  sample  with  one  of  their  g-axes 
parallel  to  the  external  magnetic  field.  Thus,  an  ENDOR  experiment  performed  at  gx  in 
figure  2-3,  will  probe  those  species  in  the  sample  with  gx  parallel  to  B0.  This  effect 
creates  a crystal-like  order  in  the  sample  at  the  edges,  gx  and  gz,  with  an  average  of  the 
various  orientations  present  between  the  extremes  of  the  powder  spectrum.  If  the 
hyperfine  coupling  is  weak,  i.e.  less  than  the  Larmor  frequency,  the  ENDOR  spectrum  is 
centered  around  the  nuclear  Larmor  frequency  of  the  coupled  nucleus,  with  ENDOR  lines 
split  with  values  equal  to  the  hyperfine  coupling  constant,  A. 


Vrf  — 


vH  ± — 
H 2 


(2.9) 


Equation  2.9  illustrates  the  relationship  between  the  frequency  of  the  external  RF,  vRF, 
and  the  proton  precessional  frequency,  vH-  This  equation  has  been  simplified  to  the 
isotropic  case,  where  the  angle  between  the  external  magnetic  field  vector  and  the  vector 
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connecting  the  electron  and  the  nucleus  is  insignificant.  So,  in  equation  2.9,  we  should 
think  of  A as  the  effective  hyperfine  coupling  constant,  rather  than  a component  of  the  A- 
tensor,  or  strictly  an  isotropic  value.  The  details  of  orientation  selection  will  be  discussed 
further  in  the  following  section.  These  experimental  details  have  been  thoroughly 
described  elsewhere  and  have  only  been  highlighted  here  [70,71], 

Orientational  Selection 

In  the  realm  of  high  fields,  the  distribution  of  molecular  orientations,  and  hence,  g- 
tensor  orientations,  are  separated  such  that  the  orientations  may  be  selected  by  changing 
the  magnetic  field  [85,86],  By  selecting  an  orientation,  a powder  or  frozen  solution, 
sample  may  be  studied  in  a partially  ordered,  crystalline,  way.  This  work  uses 
orientational  selection  to  perform  ENDOR  experiments  on  individual  orientations  to 
obtain  hyperfine  tensor  magnitudes  and  orientations  relative  to  the  g-tensor  axis  system 
[87-90].  The  success  of  this  technique,  and  its  utility  toward  organic  radicals,  has  been 
demonstrated  in  the  past  [91-93],  The  hyperfine  coupling  tensor,  A,  is  a sum  of  the 
isotropic  and  anisotropic  interaction  between  magnetic  nuclei  and  the  electron  density. 
The  isotropic  component  of  this  tensor  results  from  contact  interaction,  known  as  Fermi 
contact  interaction,  which  is  generated  through  bonds  and  is  directly  proportional  to  the 
spin  density  on  the  atom.  The  anisotropic  component  of  A results  from  dipolar 
interaction  between  the  field  produced  by  the  nucleus  and  the  field  generated  by  the 
electron  density.  Equation  2.10  describes  the  dipolar  field,  Hd,  induced  by  an  electronic 
dipole  at  any  point  in  space. 

Hd  =^^rr-E)=M,-Ad 
r 


(2.10) 
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In  this  equation,  \ie  represents  the  electron  magnetic  moment,  r is  the  distance  from  the 

electron  to  the  nucleus,  r is  a unit  vector  connecting  the  electron  and  the  nucleus,  and  E 
is  the  unit  dyadic.  Thus,  in  frequency  units,  the  hyperfine  coupling,  Ay,  can  be  expressed 
including  the  dipolar  and  isotropic  contributions. 


The  first  term  on  the  right  side  of  equation  2.1 1 is  the  dipolar  contribution,  Ad,  with  the 
second  term  obviously  equal  to  the  isotropic  contribution,  Aiso,  which  is  non-zero  on  the 
diagonal  elements  of  the  tensor.  The  three  spatial  components  of  A lead  to  the  total 
hyperfine  interaction  tensor,  shown  in  equation  2.12. 


Since  an  ENDOR  experiment  is  primarily  a nuclear  magnetic  resonance  experiment,  the 
total  magnetic  field  at  the  nucleus  must  be  known  in  order  to  predict  the  transition 
frequencies.  The  nuclear  energy  levels  result  from  the  nuclear  Zeeman  interaction,  which 
get  perturbed  by  the  hyperfine  interaction  with  the  electronic  moment  and  are  dependent 
on  the  effective  field  direction.  With  these  relationships  accounted  for,  an  expression  for 
the  frequency  of  expected  ENDOR  transitions  can  be  written  as  equation  2.13  [86]. 


Equation  2.13  gives  the  resonant  frequency  of  a nucleus  having  known  coordinates  as  a 

function  of  the  magnetic  field  direction  [hi,  fi2,  hi]  and  magnitude.  In  this  equation,  H 
represents  the  dipolar  field  induced  by  the  electron  at  any  point  in  space.  Equation  2.13 


(2.11) 


(2.12) 


(2.13) 
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shows  that,  contrary  to  equation  2.8,  the  nuclear  frequencies  are  dependent  on  the 
electronic  g-value  and  that  the  ENDOR  transitions  are  not  necessarily  symmetric  about 
the  nuclear  Larmor  frequency.  The  hyperfine  coupling  expression  can  be  simplified  by 
assuming  that  the  electron  and  nuclear  moments  are  quantized  along  the  direction  of  the 
applied  external  field,  and  that  the  spin-orbit  interactions  have  a negligible  effect  on  the 
g-value.  Equation  2.14  results  from  these  assumptions  [86]: 

A = (3  cos 2 0-1)  +Aiso  (2.14) 

r 

This  equation  relates  the  effective  hyperfine  splitting  to  the  dipolar  and  isotropic 
hyperfine  couplings  for  any  dipolar  angle,  0. 

Pulsed  EPR 

Pulsed  experiments  have  many  functions,  but  in  this  work  they  have  been  used  to 
identify  nuclei  in  the  sample.  These  experiments  hold  the  B0  field  constant  while  varying 
the  power  of  the  magnetic  component  of  the  microwave  radiation,  Bj.  Figure  2-5 
illustrates  the  effect  of  microwave  pulses  on  the  net  spin  vector,  S.  The  resultant  of  this 
vector’s  motion  is  the  bulk  magnetization,  M.  This  bulk  magnetization  can  be  composed 
of  many  different  groups  of  electron  spins,  called  spin  packets.  Each  packet  represents 
magnetic  moments  under  a slightly  different  magnetic  environment  within  the  sample. 
Figure  2-5(a)  depicts  the  motion  of  a spin  moment  in  the  presence  of  an  external 
magnetic  field.  The  resulting  vector  component  of  this  moment  is  M,  shown  in  figure  2- 
5(b).  Frame  (c)  demonstrates  the  effect  of  a 90°  microwave  pulse  on  the  net 
magnetization.  Frames  (d)  through  (f),  adopt  the  rotating  frame  point  of  reference.  Since 
the  net  spin  vector,  M,  is  precessing  about  B0,  we  can  simplify  the  picture  by  leaving  the 
laboratory  frame  and  jumping  on  the  x-  or  y-axis  of  the  precessing  frame.  Upon 
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application  of  a microwave  pulse,  Bj,  perpendicular  to  B0,  M undergoes  changes  in 
orientation  relative  to  B0.  The  actual  displacement  angle  is  determined  by  the  length  and 
the  power  of  the  pulse,  and  is  known  as  the  Rabi  angle. 


In  equation  2.15,  £2,  is  the  Rabi  angle,  t is  the  length  of  the  microwave  pulse,  and  vr  is 
the  Rabi  frequency.  The  Rabi  frequency  is  given  by: 


The  frequency  of  the  microwave  field  Bi  is  given  by  Vi,  and  at  resonance,  v = Vb.  Once 
displaced  through  £2,  the  spin  packets  will  dephase  as  a result  of  inhomogeneities  in  the 
effective  magnetic  field  (figure  2-5(d)).  This  dephasing  time  is  called  T2  or  the 
transverse  relaxation  time  and  is  inversely  related  to  the  strength  of  the  spin-spin 
coupling.  This  dephasing,  or  relaxation,  can  result  from  homogeneous  or  inhomogeneous 
effects.  Homogeneous  relaxation  occurs  when  the  magnetic  field  is  the  same  throughout 
the  excited  spins.  Inhomogeneous  relaxation  occurs  when  the  magnetic  field  felt  by  the 
spins  varies  throughout  the  sample.  This  inhomogeneity  can  result  from  an 
inhomogeneous  external  field,  unresolved  hyperfine  structure,  anisotropic  hyperfine  and 
Zeeman  interactions,  and  dipolar  interactions  [69],  These  effects  alter  the  effective 
magnetic  field  felt  by  the  individual  spins,  affecting  their  interaction  with  the  external 
radiation.  This  dephasing  can  be  reversed  by  refocusing  the  spin  packets  with  a 180° 
pulse.  The  refocusing  of  the  spin  packets  generates  a spin  echo  (Figure  2-5(e)).  This 
effect  was  first  discovered  by  Hahn,  who  used  a 90°  pulse,  followed  by  a 180°  pulse  [85]. 
Hahn  echos  are  created  with  two  pulses,  but  echos  can  be  induced  with  many  pulses,  as 
long  as  the  magnetization  gets  refocused  in  the  xy  plane.  As  long  as  the  net 


£2  = 27i(x)(vr) 


(2.15) 


(2.16) 
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Figure  2-5  Effects  of  microwave  pulses  on  (a)  the  spin  magnetic  moment,  S,  and  (b-f)  the 
net  magnetization,  M.  The  spins  are  depicted  in  the  laboratory  (stationary) 
frame  in  (a)  through  (c).  (d)  through  (f)  assume  the  rotating  frame  of 
reference  where  point  of  view  follows  the  precession  of  the  magnetization. 

magnetization,  M,  lies  in  the  xy  plane,  the  echo  will  be  detected.  The  spin  packets  will 
relax  to  their  equilibrium  orientation,  0 = 0,  with  time  constant,  Tj.  Different  pulse 


sequences  can  be  used  to  determine  T2  and  Tj,  which  are  used  to  optimize  the  parameters 
for  pulsed  ENDOR  and  ESEEM  experiments. 

Pulsed  ENDOR 

Pulsed  ENDOR  experiments  allow  the  detection  of  NMR  through  the  intensity  of 
the  electron  spin  echo.  One  advantage  over  cw-ENDOR  is  the  detection  of  the  ENDOR 
effect  in  the  absence  of  all  incident  radiation.  That  is,  the  microwave  and  RF  fields  are 


turned  off  during  the  detection.  Pulsed  ENDOR  experiments  manipulate  the  spin  packets 
with  microwave  pulses,  then  measure  the  affect  of  an  RF  pulse  on  the  bulk 
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magnetization.  Davies  [94]  and  Mims  [95]  developed  two  variations  on  the  pulsed 
ENDOR  experiment. 

Davies  ENDOR  applies  a selective,  saturating  microwave  pulse,  stimulates  NMR 
transitions  with  an  RF  pulse,  then  refocuses  the  resulting  magnetization  with  a Hahn  echo 
sequence.  The  pulse  sequence  for  Davies  ENDOR  is  shown  in  figure  2-6.  The  first  180° 
pulse  fully  inverts  the  spin  population  of  one  of  the  EPR  transitions,  such  that  the  bulk 
magnetization  is  inverted  along  B0.  This  pulse  saturates  a portion  of  the 
inhomogeneously  broadened  EPR  line.  This  “hole”  in  the  EPR  line  is  then  filled  in  as 
NMR  transitions  are  induced.  The  NMR  transition  moves  spin  population  back  into 


RF 


Figure  2-6  Davies  ENDOR  pulse  sequence.  A saturating  microwave  7t  pulse  is  followed 
by  a 7t  RF  pulse,  then  a Hahn  echo  sequence,  7t/2  - n,  refocuses  the 
magnetization.  The  intensity  of  the  spin  echo  generated  will  change  if  nuclear 
transitions  are  induced  by  the  RF  pulse. 

the  hole  created  by  the  saturating  microwave  pulse.  So,  Davies  ENDOR  detects  the 
restoration  of  spin  echo  intensity  as  a function  of  RF  frequency.  The  bandwidth  of  the 
initial  microwave  pulse  should  be  much  less  than  the  hyperfine  coupling  constant,  as 
shown  in  equation  2.17. 


A 
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v=r<<a‘ 

fJW 


(2.17) 


In  this  equation,  Av  represents  the  bandwidth  of  the  initial  microwave  pulse,  tHW 


represents  the  duration  of  the  initial  microwave  pulse,  and  aiS0  is  the  isotropic  hyperfine 
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coupling  constant.  If  the  range  of  frequencies  in  the  pulse  is  too  large,  on  the  order  of 
aiS0,  then  the  range  of  spin  packets  affected  by  the  microwave  pulse  may  overlap  with  the 
hyperfine  levels,  diminishing  the  ENDOR  effect,  since  the  levels  that  fill  in  the  hole  will 
be  saturated  as  well. 

Mims  ENDOR  applies  a series  of  three  non-selective  90°  pulses  with  the  second 
and  third  microwave  pulses  separated  by  an  RF  pulse  (see  figure  2-7).  The  first  two 
microwave  pulses  create  a grated  polarization  pattern  between  the  electron  spins,  Sz,  and 
the  two-spin  order  terms,  SZIZ.  If  an  RF  pulse  induces  NMR  transitions,  the  electron  spins 
coupled  to  those  nuclear  spins  (SZIZ)  will  be  dephased  from  the  bulk  magnetization,  thus 
reducing  the  intensity  of  the  electron  spin  echo.  This  technique  is  more  sensitive  to  small 
hyperfine  coupling  constants  (a  < 4 MHz)  but  is  also  known  to  have  “blind  spots.”  These 
blind  spots  result  from  incomplete  refocusing  of  the  bulk  magnetization.  The  application 
of  a 7t/2  microwave  pulse  along  the  y-axis  of  the  rotating  frame  affects  the  bulk 
magnetization  along  the  x-axis.  The  remaining,  dephased  spin  packets  in  the  xy-plane 
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Figure  2-7  Mims  ENDOR  pulse  sequence.  Two  n!2  microwave  pulses  are  applied, 

followed  by  a n RF  pulse,  then  another  nJ2  microwave  pulse  to  refocus  the 
spins.  If  the  RF  pulse  “flips”  any  nuclear  spins  that  are  coupled  to  the  electron 
spin,  these  coupled  spin  packets  will  not  contribute  to  the  spin  echo. 

are  rotated  90°  about  the  y-axis,  since  the  microwave  pulses  are  applied  along  the  y- 

direction.  So,  the  spin  packets  on  the  y-axis  are  not  refocused  in  the  echo,  which  reduces 
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the  echo’s  intensity.  The  relationship  between  the  pulse  separation  and  the  hyperfine 
interaction  is  shown  in  equation  2.18. 


t = 


2m 


(2.18) 


In  this  equation,  x represents  the  time  between  the  first  and  second  microwave  pulses,  n is 
an  integer,  and  ajS0  is  the  isotropic  hyperfine  coupling  constant  given  in  units  of  27tHz. 
Because  of  the  condition  shown  in  equation  2.18,  Mims  ENDOR  data  should  be  taken  at 
various  x values. 


ESEEM 

Electron  Spin  Echo  Envelope  Modulation  (ESEEM)  was  first  demonstrated  by 
Rowan  et  al.  in  1965  [96].  Mims  also  contributed  to  the  theory  and  experimental 
techniques  for  this  experiment  in  1975  [97].  ESEEM  is  a valuable  technique  used  to 
locate  interacting  nuclei  [101]  and  to  measure  hyperfine  couplings  of  distant  nuclei  or 
nuclei  with  small  nuclear  moments.  This  experiment  is  essentially  a spin  echo 
experiment  with  a varied  pulse  separation.  This  technique  capitalizes  on  the  relative 
strength  of  the  electron  and  nuclear  magnetic  moments.  Since  the  nuclear  moment  is 
approximately  three  orders  of  magnitude  smaller  than  the  electron’s,  the  response  of 
these  magnetic  moments  to  microwave  pulses  varies  by  a few  orders  of  magnitude  in  the 
time  domain.  That  is,  the  microwave  pulse  reorients  the  electron’s  moment  much  more 
quickly  than  the  nuclear  moment.  The  effect  of  this  lag  time  between  the  moments 
results  in  incomplete  refocusing  of  the  spin  packets,  which  varies  the  intensity  of  the 
echo.  If  the  microwave  pulse  interval,  T,  is  divisible  by  the  precessional  periods  of 
coupled  nuclei,  the  intensity  of  the  spin  echo  will  be  enhanced.  The  ESEEM  effect  can 
be  observed  with  two  microwave  pulses,  but  stimulated  (three  pulse)  ESEEM  provides 
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greater  frequency  resolution  and  eliminates  hyperfine  artifacts  that  result  from  phase 
interference  between  nuclei  [98].  The  ESEEM  spectrum  is  a plot  of  T versus  the 
intensity  of  the  spin  echo,  where  T is  the  varied  pulse  separation  (see  figure  2-8).  As  T is 
increased,  nuclear  couplings  modulate  the  first  order  decay  of  the  spin  echo,  which  has 
time  constant  T2  for  2-pulse  ESEEM  and  Tj  for  3-pulse  ESEEM.  The  ESEEM  effect 
arises  from  the  magnetic  fields  generated  by  coupled  magnetic  nuclei,  which  give  rise  to 
variation  in  the  precessional  frequency  of  the  spin  packets.  Since  the  effect  of  the  nuclear 
field  on  the  electron  is  a function  of  the  spin  state  of  the  electron,  the  spin  packets 
coupled  to  nuclei  can  have  a slightly  greater  precessional  frequency  (if  the  nuclear  field  is 
additive)  or  slightly  smaller  precessional  frequency  (if  the  nuclear  field  is  subtractive) 


(b) 


Figure  2-8  (a)  2 pulse  and  (b)  3 pulse  ESEEM  sequences.  The  delay  time,  T,  is  increased 
throughout  the  experiment  and  the  response  of  the  echo  is  detected  at  each  T 
value. 

than  the  bulk  magnetization.  This  fragments  the  magnetization  and  results  in  a phase 
difference  between  the  coupled  packets  and  the  non-coupled  packets.  The  interference 
between  the  precession  of  the  coupled  packets  with  the  bulk  magnetization  results  in  a 
modulation  of  the  echo  intensity,  which  is  ESEEM. 

The  ESEEM  effect  can  also  be  explained  through  the  allowance  of  forbidden 
transitions.  In  order  to  fully  explain  this  phenomenon,  the  classical  picture  must  be 
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enhanced  with  quantum  effects.  Since  the  microwave  pulse  covers  a range  of  frequencies 
that  encompass  the  Mi  levels  of  a given  spin  sublevel,  the  wavefunctions  become  mixed. 
At  small  hyperfine  couplings,  formally  forbidden  transitions  become  weakly  allowed, 
such  as  (AMs  = ±1,  AMi  = + 1)  and  can  interfere  with  the  allowed  (AMs  = ±1,  AMi  = + 0) 
transitions.  This  interference,  or  branching,  modulates  the  echo  intensity  with  a 
frequency  dependent  on  the  spacing  of  the  nuclear,  Mi,  sublevels  [99,100].  The 
probability  of  this  branching  is  determined  by  the  forbidden  transitions  that  become 
allowed  if  the  non-selective  microwave  pulse  flips  the  nuclear  spin,  which  creates  a two 
spin  order.  If  branching  does  not  occur,  the  upper  and  lower  states  become  coherent  and 
generate  a spin  echo.  The  branching  then  introduces  coherent  relationships  between  the 
upper  and  lower  manifolds,  which  modulates  the  echo  intensity. 


CHAPTER  3 

EXPERIMENTAL  PROCEDURES 

Bacteriochlorophyll  a Purification 

Crude  BChl  a was  harvested  from  the  cells  of  Rhodobacter  sphaeroides  [102],  The 
contents  of  the  bacterial  cells  were  separated  using  column  chromatography.  A 
powdered  sugar  column  was  used  for  the  purification  [67,103].  A 41mm  Kontes  column 
(K420530-0273)  was  packed  with  16”  of  powdered  sugar  (1”  ID).  The  sample  (100  mg) 
was  loaded  onto  the  column  in  a minimal  amount  of  10%  acetone  in  1:9  ethyl 
ether/hexane  (v/v).  Elution  was  initiated  with  200  mL  aliquots  of:  1:9  ethyl 
ether/hexane,  2%  acetone  in  1:9  ethyl  ether/hexane,  3%  acetone  in  1:9  ethyl  ether/hexane, 
5%  acetone  in  1:9  ethyl  ether/hexane,  and  finally  7%  acetone  in  1:9  ethyl  ether/hexane. 
The  eluted  species  were  monitored  using  thin  layer  chromatography  on  silica  plates 
developed  in  70:30  hexane/acetone.  The  column  was  run  in  the  dark,  or  in  the  presence 
of  indirect  green  light.  The  BChl  a fractions  were  combined,  evaporated  to  dryness  using 
a Buchi  Rotovapor  R-200  pumped  with  a Buchi  Vac  V-500,  and  stored  under  argon,  in 
the  dark,  at  -20°C.  The  fractions  were  also  monitored  with  UV-Visible  Spectroscopy  on 
a Varian  Cary  500e  Spectrometer.  Figure  3-1  contains  the  UV-Vis  spectra,  which  are  in 
agreement  with  the  expected,  published  data. 


32 


33 


Figure  3-1  UV-Vis  spectra  of  purified  BChl  a in  acetone  and  1:9  ethyl  ether/hexane. 

Bacteriopheophytin  a Generation 

BPheo  a was  generated  from  bacteriochlorophyll  a.  BChl  a was  purified  as 
described  in  the  previous  section.  Once  pure,  BPheo  a was  dissolved  in  glacial  acetic 
acid,  which  protonates  the  two  nitrogens  coordinated  to  the  magnesium  atom,  hence 
freeing  the  magnesium  ion  from  the  molecule.  The  bacteriopheophytin  was  extracted 
using  a phase  separation  technique  with  toluene  and  water.  The  bacteriopheophytin  goes 
into  the  toluene  layer,  while  the  magnesium  ions  and  acetic  acid  are  in  the  aqueous  layer. 
The  toluene  layer  was  repeatedly  washed  with  distilled  water,  and  removed  from  the 
aqueous  portions.  The  toluene  layer  was  then  evaporated  under  reduced  pressure  and  the 
product  was  examined  with  a Cary  500e  UV-Visible  Spectrometer,  spectrum  shown  in 
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figure  3-2.  An  alternative  procedure  developed  by  Lotjonen  and  Hynninen  in  1983  [103] 
can  also  be  used  for  demetallation. 


Figure  3-2  UV-Vis  absorption  spectrum  of  BPheo  a in  acetone  and  1:9  ethyl 
ether/hexane. 

Chlorophylls  a and  b — Extraction  and  Purification 

Plant  chlorophylls  a and  b were  isolated  from  spinach  leaves  using  a modified 
dioxane  method  [104,105].  In  the  presence  of  dim  green  light,  60g  of  spinach  leaves 
were  ground  up  in  a mortar  and  pestle  in  the  presence  of  90mL  acetone  and  0.30g 
NaHP04  7H2O.  This  leaf  mixture  was  filtered  through  a coarse  Buchner  funnel  into  a 
300mL  round  bottom  flask.  60mL  dioxane  was  added,  then  distilled  water  was  added 
dropwise  until  precipitation  became  visible  (~30mL).  The  flask  was  placed  in  the  freezer 
(-20°C)  for  45  minutes.  The  fully  precipitated  mixture  was  then  centrifuged  at  15,000 
relative  centrifugal  force  (RCF)  for  30  minutes  at  5°C  using  an  Eppendorf  5804R 
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centrifuge.  The  pellet  was  dissolved  in  ethanol  and  evaporated  to  dryness  under  reduced 
pressure.  This  pellet  contains  carotenoids,  Chi  a and  Chi  b.  These  components  were 
separated  with  column  chromatography  using  a custom  made  column;  48”  tall,  2”  ID. 
This  column  was  packed  with  40”  of  powdered  sugar  for  80mg  of  crude  leaf  extract.  The 
sample  was  loaded  in  10%  acetone  in  1:9  ethyl  ether/hexane.  The  carotenoids  were 
eluted  with  hexane  and  1:9  ethyl  ether/hexane.  A solvent  gradient  was  established  with 
200  mL  aliquots,  from  1%  acetone  in  1:9  ethyl  ether/hexane  to  10%  acetone  in  1:9  ethyl 
ether/hexane.  Chi  a was  eluted  with  3%  acetone  in  1:9  ethyl  ether/hexane,  and  Chi  b was 
eluted  with  5%  acetone  in  1:9  ethyl  ether/hexane.  This  column  was  run  in  total  darkness 
and  the  fractions  were  monitored  with  TLC  and  with  a Cary  500e  UV-Visible 


Figure  3-3  Chlorophyll  a (solid)  and  b (dotted)  UV-Vis  absorption  spectra  in  1:9 
ether/hexane. 
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spectrometer  (Varian  Instruments).  The  UV-Vis  spectra  are  shown  in  figure  3-3  for  both 
pigments  in  1:9  ether/hexane.  The  pure  components  were  stored  under  argon,  in  the  dark, 
at  -20°  C. 

Radical  Anion  Generation 

Approximately  1-2  mg  of  the  pigment  was  dissolved  in  acetone  (C3H6O)  or 
methylene  chloride  (CH2CI2)  and  added  to  the  reaction  vessel  shown  in  figure  3-3,  which 
had  been  dried  in  a 150°C  oven  overnight.  The  solvent  was  then  evaporated  under 
reduced  pressure,  ~104  torr.  Potassium  was  introduced  into  the  vessel  in  glass  capillaries, 
which  were  placed  in  the  potassium  arm.  After  repeated  flame  drying  of  the  entire  vessel 

Vacuum 


t 


Figure  3-3  Anaerobic  reaction  vessel  made  from  10mm  borosilicate  tubing  with  4mm  or 
5mm  quartz  EPR  tubing  and  a borosilicate-quartz  graded  seal.  The  tapered 
sections  of  the  vessel  facilitate  glassblowing  under  high  vacuum. 
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(with  care  taken  to  keep  the  pigment  in  the  sample  arm  around  30°C),  the  potassium  arm 
of  the  reaction  vessel  was  heated  to  produce  a potassium  mirror.  This  mirror  was  made 
several  times  over  several  hours,  to  allow  sufficient  reaction  with  any  residual  water  in 
the  system.  The  tip  of  the  potassium  arm  was  eventually  sealed  off  with  a torch  to  isolate 
any  byproducts  from  the  melting  and  boiling  of  the  potassium.  The  solvent  was  then 
distilled  over  from  another  flask,  which  contained  the  freshly  distilled  and  degassed 
solvent  over  a sodium  mirror.  This  solvent  was  attracted  to  the  reaction  vessel  by 
submerging  the  sample  arm  in  liquid  nitrogen.  The  reaction  vessel  was  then  sealed  off 
from  the  vacuum  system  and  once  cool,  the  reduction  was  achieved  by  pouring  the 
solution  over  the  mirror.  This  reaction  was  monitored  with  EPR,  and  once  a sufficient 
concentration  of  the  product  was  formed,  the  EPR  arm  and  the  sample  arm  were  frozen 
simultaneously  and  each  sealed  off  with  an  O2/H2  mixed  gas  torch. 

X-Band  EPR/ENDOR 

The  X-band  EPR  was  done  on  a Bruker  Elexsys  e580  pulsed  spectrometer.  This 
instrument  uses  an  e580-1010  pulsed  microwave  bridge  with  a Gunn  microwave  source 
and  an  Applied  Systems  Engineering  117x  TWT  amplifier  to  generate  high  powered 
pulses.  The  static  magnetic  field  is  maintained  by  a 1.5  T electromagnet.  The  data  was 
collected  and  stored  on  a Silicon  Graphics  O2  Workstation,  using  the  supplied  Bruker 
Xepr  software.  An  ENI  A-300  Radio  Frequency  amplifier  was  used  to  generate  and 
amplify  the  RF  pulses  required  for  pulsed  ENDOR  experiments.  The  cw  experiments 
used  a TE101  rectangular  cavity  with  an  Oxford  CF-900  continuous  flow  cryostat.  The 
pulsed  experiments  were  performed  using  a Flexline  Bruker  ENDOR  cavity  with  an 
Oxford  CF-935  helium  flow  cryostat.  Cryogenic  temperatures  were  maintained  by  an 
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Oxford  ITC  503  Temperature  Control  Unit,  an  Oxford  GFS-650  Transfer  Tube,  and  an 
Oxford  VC40  Gas  Flow  Controller. 

Transmission  High  Field  cw-EPR 

The  transmission  high  field  cw-EPR  experiments  were  performed  on  a home-built 
setup  at  NHMFL  that  has  been  previously  described  [98].  This  “single  pass”  setup  is 
depicted  in  figure  3-4.  Microwaves  are  generated  using  a 93  GHz  Gunn  source,  which 
was  coupled  to  a Shottky  diode  harmonic  generator  and  high-pass  filters,  which  allows 
operation  at  frequencies  up  to  550  GFIz.  The  Gunn  source  is  coupled  to  a YIG  transistor 
oscillator  to  lock  the  frequency.  This  microwave  radiation  is  sent  through  oversized, 
cylindrical,  brass  waveguides  to  the  sample  and  then  to  a liquid  helium  cooled,  QMC 
In:Sb  hot  electron  bolometer  detector  (Model  QFI/3BI).  The  absorption  of  microwaves  is 
measured  directly,  without  a reference  arm  or  a cavity.  The  sample  probe  contains  a 
modulation  coil  that  was  operated  at  1 kHz,  with  modulation  amplitudes  between  1 Gauss 
and  10  Gauss.  This  instrument  is  outfitted  with  an  Oxford  Instruments  Teslatron  17 
Tesla,  superconducting  magnet,  which  uses  a sweep  coil  for  sweeps  less  than  0.1  Tesla. 
The  EPR  signal  is  monitored  with  a millimeter  Vector  Network  Analyzer  (MVNA), 
Model  MVNA-8-350,  for  signal  processing.  The  field  modulation  was  generated  using  a 
Thurlby  & Thandar  function  generator  (Model  TG210)  which  was  amplified  with  a 
MOS-FET  current  amplifier  (Model  X-820HD).  The  modulation  coil  covers  a 50  mm 
region  of  the  waveguide  and  is  tightly  wound  about  a G-10  mandrel.  Lock-In  amplifiers 
from  Stanford  Research  Systems,  Inc.  (Model  DSP  SR830)  were  also  employed  to  isolate 
the  signal  component  that  follows  the  modulation  frequency.  Cryogenic  temperatures 
were  achieved  through  the  use  of  an  Oxford  Instruments  CF-1200  continuous  flow 


cryostat. 
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Microwave 
Sources 
(95-330  GHz) 


Main  Coil:  0-1 7 Tesla 
Sweep  Coil:  +/-  0. 1 Tesla 

Figure  3-4  Transmission  cw-EPR  block  diagram.  Figure  adapted  from  [98], 

Quasioptic  High  Field  cw-EPR/ENDOR 

The  quasi-optic  high  field  EPR  and  ENDOR  experiments  were  done  on  a home- 

built  spectrometer  at  NHMFL  [106],  This  is  a superheterodyne  instrument,  shown  in 

figure  3-5,  uses  two  Gunn  sources  (G1  and  G2)  and  two  Shottky  diode  detectors  (Ml  and 

M2),  mixed  with  a quasi-optical  bridge  (dotted  region).  As  the  microwaves  (from  Gl) 

enter  the  bridge,  they  are  attenuated  with  a wire-grid  polarizer,  and  a Faraday  rotator. 
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The  quasi-optical  bridge  is  composed  of  elliptical  aluminum  mirrors  that  mix  and  guide 
the  microwaves  to  a Martin-Pupplet  interferometer,  which  converts  the  polarization 
before  entering  the  oversized  corrugated,  cylindrical  waveguide  in  the  sample  probe.  The 


SI  S2 


Figure  3-5  Quasi-optical  superheterodye  instrumental  setup  [106], 
microwaves  return  from  the  sample,  get  split  into  two  components  by  a wire-grid 
polarizer.  The  co-polarized  beam  returns  to  the  circulator  and  then  to  the  reference 
Schottky  diode  (Ml),  while  the  cross-polarized  component  goes  to  the  signal  detector 
(M2).  The  returning  microwaves  are  mixed  with  the  second  Gunn  source  (G2),  operating 
6 GHz  below  the  main  Gunn,  and  provides  a local  oscillator  signal  for  the  Shottky 
diodes.  At  the  sample,  a corrugated  horn  is  used  to  focus  the  microwaves  down  to  the 
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sample  holder,  which  is  made  of  Rexalite,  shown  in  figure  3-6.  This  sample  holder 
accommodates  5 mm  quartz  tubes  with  length,  less  than  60  mm.  This  setup  allows  for 


Figure  3-6  Quasi-optical  sample  holder  for  frozen  liquid  samples,  for  240  GHz  cw- 
EPR/ENDOR  measurements. 

cold  loading  of  the  sample  into  the  magnet.  A modulation  coil  is  fitted  beneath  the 
holder,  and  a 6 turn  ENDOR  coil  is  oriented  around  the  sample  tube.  A Rhode  and 
Schwartz  RF  source  generates  the  radiofrequency  radiation,  which  is  amplified  by  an 
Amplified  Research  RF  amplifier. 

Magnetic  Field  Calibration  Procedure 

The  uncertainty  in  the  magnetic  field  was  reduced  through  the  use  of  standards  with 
well  known  Zeeman  and  hyperfine  parameters.  The  two  standards  employed  in  this  work 
were  MgO  containing  0.05%  Mn  ( I = 5/2  ) and  silicone  doped  with  phosphorus  (1  = 1/2  ). 
One  of  these  standards  was  included  with  the  sample  in  the  instrument  and  the  position  of 
the  standard’s  lines  provide  a guide  to  correct  the  field.  The  gross  uncertainty  in  the  field 
arises  from  the  hysteresis  in  the  current  within  the  sweep  coil  used  to  sweep  the  field 


during  the  EPR  measurements. 
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The  manganese  standard,  MgO  powder,  was  incorporated  into  epoxy,  allowed  to 
harden,  then  cut  to  fit  into  the  waveguide.  The  EPR  spectrum  was  measured  such  that  the 
Mn  lines  were  resolved  and  intense  enough  to  simulate.  The  field  was  corrected  by 
fitting  the  manganese  lines  in  the  experimental  data  using  EPR  simulation  software, 
taking  the  ratio  of  the  experimental  g-value  to  the  accepted  g-value,  then  shifting  the 
spectrum  by  this  correction  factor.  The  experimental  Mn  spectrum  was  fitted,  and  the  g- 
value  for  this  fitting,  g(Mn)exp,  was  then  divided  by  the  accepted  value,  g(Mn)acc  = 
2.00101  + 5xl0"5.  The  ratio  of  g(Mn)exp/g(Mn)acc  gives  the  correction  factor.  Once  the 
powder  pattern  was  fitted,  using  the  DOUBLET  software,  the  experimental  g-values  for 
the  sample  were  divided  by  the  correction  factor.  This  procedure  shifts  the  g-values, 
which  are  linearly  proportional  to  the  magnetic  field  values. 

For  the  bacteriopheophytin  a radical  anion  samples,  the  field  calibration  procedure 
was  performed  using  a Si:P  standard.  This  procedure  is  identical  to  the  Mn  field 
correction.  The  experimental  phosphorus  g-value,  g(P)exp,  was  determined  by  fitting  the 
experimental  phosphorus  lines,  which  was  then  compared  to  the  accepted  value,  g(P)aCc= 
1.998525  + 5x10 5.  The  ratio  of  these  g-values,  g(P)exp/g(P)acc,  gives  the  correction  factor 
by  which  the  experimental  powder  pattern  g-values  were  divided.  The  simulation 
procedure,  which  accounts  for  second  order  effects,  accurately  fit  the  manganese 
spectrum,  which  gives  an  experimental  g-value  for  Mn  as  well  as  indicating  linearity  in 
the  magnetic  field  throughout  the  spectrum. 

The  magnetic  field  on  the  Bruker  e580  X-band  spectrometer  was  calibrated  using 
the  perylene  radical  cation,  which  was  generated  in  concentrated  sulfuric  acid.  The 
accepted  g-value  for  this  radical  is  known  to  be  g(perylene+)iso  = 2.002569  ±6xl0'6  [69], 
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The  experimental  spectrum  was  fitted  using  the  DOUBLET  software,  and  the  ratio  of  the 
experimental  g-value  and  the  accepted  g-value  was  used  as  a correction  factor  to  shift  the 
X-band  spectra.  The  experimental  g-value  was  measured  to  be  g(perylene+)  exp  = 
2.003877  + lxlO'6,  which  corresponds  to  a shift  of  -2.20  G. 


CHAPTER  4 
RESULTS 

Bacteriochlorophyll  a Radical  Anion 

X-Band  (cw  and  pulsed)  and  high  field  (cw)  measurements  were  used  to 
investigate  the  bacteriochlorophyll  a radical  anion.  The  X-Band  data  provide  only 
isotropic  g and  A-values,  which  serve  as  a reference  for  experiments  at  higher  fields  since 
g and  A-tensor  data  are  sought  out.  The  magnitude  of  the  hyperfine  coupling  dictates  the 
linewidth  at  X and  W-bands,  but  full  resolution  in  g is  achieved  around  240  GElz.  The  X- 
band  EPR  experiments  provide  partially  resolved  isotropic  hyperfine  couplings.  Figure 
4-1  includes  experimental  (solid  line)  and  simulated  (dashed  line)  EPR  spectra.  This  X- 
band  data  shows  the  presence  of  hyperfine  coupling  interactions  that  were  estimated 
using  the  DOUBLET  program,  kindly  provided  by  Dr.  Andrew  Ozarowski.  The 
simulation  in  figure  4-1  was  generated  by  considering  contributions  from  the  three 
methine  bridge  protons  (positions  a,  (3,  5),  six  methyl  protons  (three  from  the  methyl  on 
carbon  1,  and  three  from  the  methyl  at  carbon  5),  and  two  nitrogens  (positions  22  and 
24).  Figure  1-3  presents  the  labeling  assignments  for  these  atoms.  The  hyperfine 
couplings  are  as  follows:  Amethyi  = 3.0  G (8.4  MHz),  Amethine  = 2.4  G (6.7  MHz),  and  An 
= 2.2  G (6.2  MHz).  The  EPR  line  is  centered  at  gjS0  = 2.003250  + 6xl0~6  with  a linewidth 
of  1.2  G.  This  experiment  gives  partially  resolved,  isotropic  hyperfine  coupling  constants 
since  the  liquid  state  averages  out  the  anisotropic  contribution.  This  is  in  agreement  with 
previously  published  results  [56,107],  Cryogenic  X-Band  EPR  yields  a single, 
unresolved  line  with  a linewidth  of  18  G. 
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Figure  4-1  X-Band  EPR  of  BChl  a anion  taken  at  298  K and  9.405913  GHz,  with  1.0 
mW  microwave  power,  0.5  G modulation  amplitude  at  1 kHz,  10  ms  time 
constant,  and  20  ms  conversion  time.  The  experimental  data  (solid)  was 
simulated  (dotted)  assuming  the  following  hyperfine  contributions:  6 x Amethyi 
= 3.0  G,  3 x Amethine  = 2.4  G,  and  2 x AN  = 2.2  G. 

Pulsed  experiments  were  done  at  X-Band  to  further  investigate  the  hyperfine 
interactions  in  this  species.  Mims  ENDOR  was  done  on  BChl  a radical  anion  and  shows 
two  resolved  proton  couplings  at  1.12  MHz  and  9.60  MHz  (see  figure  4-2).  Mims 
ENDOR  is  known  to  have  blind  spots  that  could  explain  the  missing  ENDOR  transitions 
expected  in  this  spectrum.  However,  these  values  fall  within  the  range  of  isotropic 
couplings  determined  by  the  X-band  cw-EPR,  shown  in  figure  4-1.  These  pulsed 
experiments  were  conducted  on  the  Bruker  e580  pulsed  EPR  spectrometer  described  in 


chapter  3. 
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Figure  4-2  X-Band  Mims  ENDOR  on  BChl  a radical  anion,  taken  at  23K  with  two  16  ns 
(90°)  pulses  separated  by  124  ns,  a 7 ps  RF  pulse,  and  a 16  ns  refocusing 
microwave  pulse. 

The  high  field  EPR  data  show  a fully  resolved  g-tensor  at  210  GHz.  Figure  4-3 
shows  a transmission  cw-EPR  experiment  done  at  210.374  GHz,  10  K,  1 mW  microwave 
power,  and  4 G field  modulation  at  1 kHz.  The  experimental  setup  for  the  transmission 
instrument  is  described  in  chapter  3.  In  addition  to  the  powder  pattern,  this  spectrum 
shows  manganese  lines  from  the  field  calibration  standard.  Three  of  the  Mn  lines  are 
visible  on  the  high  field  side  of  the  radical  signal,  one  on  the  low  field  side,  and  two 
overlap  with  the  radical  signal.  The  corrected  g-values  for  BChl  a radical  anion:  gx  = 
2.004257  ± 5xl0'5,  gy  = 2.003150  ± 5xl0'5,  and  gz  = 2.002513  ± 7xl0'5.  These  g-values 
average  to  the  isotropic  g-value,  gave  = 2.003306  + 7xl0"5. 
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Figure  4-3  (a)  High  Field  EPR  spectrum  of  BChl  a radical  anion  taken  at  10  K and 

210.374  GHz  with  ~ 0.4  mW  microwave  power  and  4 G field  modulation  at 
~1  kHz.  This  data  was  taken  on  the  transmission  setup  with  an  MgO  field 
standard  containing  manganese.  The  simulation  employed  anisotropic 
linewidths  given  by:  Ax  = 4 G,  Ay  = 5 G,  and  Az  = 1 1 G.  (b)  Temperature 
dependence  EPR  study  of  BChl  a radical  anion  at  326.616  GHz  with  ~ 0.5 
mW  microwave  power  and  5 G field  modulation  at  1 kHz. 
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High  field  EPR  and  ENDOR  measurements  were  taken  at  9 T / 240  GHz,  on  the 
superheterodyne  setup.  This  instrumental  setup  is  discussed  in  chapter  3 [106],  The 
saturated  EPR  signal  is  shown  in  figure  4-4,  with  a fully  resolved  rhombic  g-tensor. 
Figure  4-4  is  a plot  of  the  EPR  intensity  as  the  magnetic  field  is  swept  with  a sweep  coil, 
with  the  main  coil  in  persistent  mode  at  8.560  T.  This  spectrum  appears  as  an  absorption 
due  to  rapid  passage  effects  and  lifetime  broadening  which  results  from  saturation.  'H- 
ENDOR  spectra  were  taken  in  10  G steps  across  the  EPR  envelope.  ENDOR  spectra  for 
the  principle  g-values  are  shown  in  figures  4-5  through  4-7.  The  grid  lines  on  figure  4-4 
indicate  the  magnetic  field  value  for  each  ENDOR  spectrum.  The  principle  g-values  fall 
around  -7  mT  (gx),  -1.5  mT  (gy),  and  +2  mT  (gz),  which  correspond  to  8.5988  T (gx), 
8.6036  T (gy),  and  8.6067  T (gz)  at  239.96  GHz.  The  ENDOR  spectra  have  been 


Field  Offset  (mT) 
Relative  to  8.560T 


Figure  4-4  Saturated  BChl  a anion  EPR  signal  with  fully  resolved  g-tensor.  The 

spectrum  was  taken  at  30  K and  239.96  GHz,  with  8.560  T persistent  field, 
and  5 G field  modulation  at  1kHz. 
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simulated  using  GENDOR  [107]  to  estimate  the  number  and  magnitude  of  each  pair  of 
ENDOR  transitions.  Figure  4-9  shows  a stack  plot  of  all  the  BChl  a radical  anion 
ENDOR  spectra.  These  spectra  change  as  a function  of  g,  indicating  anisotropic 
behavior,  and  hence,  orientational  selection.  The  effective  hyperfine  splittings  obtained 
from  the  simulations  are  shown  in  figure  4-9.  These  data  show  slight  anisotropy  in  the 
hyperfine  couplings  of  certain  protons.  The  0.20  MHz  splitting  seems  to  be  invariant 
throughout  the  ENDOR  spectra.  This  splitting  may  be  a natural  void  in  the  spectra 
resulting  from  the  simultaneous  excitation  of  both  hyperfine  transitions.  This  effect 
prevents  nuclear  polarization,  which  must  be  achieved  to  see  an  absorption  of  energy.  In 
addition,  the  RF  frequency  is  modulated  at  40  kHz  during  the  ENDOR  experiments.  This 
modulation  attempts  to  excite  more  spins,  but  since  the  detection  does  not  follow  this 
modulation,  the  experiment  is  blind  to  features  less  than  80  kHz  wide. 


VRF  - VH  <MHZ) 

Figure  4-5  High  field  cw-ENDOR  spectrum  of  BChl  a radical  anion,  measured  at  30K, 
with  g = 2.00402  (-7  mT).  The  experimental  spectrum  (solid)was  simulated 
(dotted)  using  the  GENDOR  software  to  estimate  the  hyperfine  tensors.  The 
effective  hyperfine  splittings  are  plotted  in  figure  4-9. 
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Figure  4-6  High  field  cw-ENDOR  spectrum  of  BChl  a radical  anion,  measured  at  30K, 
with  g = 2.00303  (-2  mT).  The  experimental  spectrum  (solid)was  simulated 
(dotted)  using  the  GENDOR  software  to  estimate  the  hyperfine  tensors.  The 
effective  hyperfine  splittings  are  plotted  in  figure  4-9. 
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Figure  4-7  High  field  cw-ENDOR  spectrum  of  BChl  a radical  anion,  measured  at  30K, 
with  g = 2.00218  (+2  mT).  The  experimental  spectrum  (solid)was  simulated 
(dotted)  using  the  GENDOR  software  to  estimate  the  hyperfine  tensors.  The 
effective  hyperfine  splittings  are  plotted  in  figure  4-9. 
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Figure  4-8  Stack  plot  of  BChl  a radical  anion  ENDOR  spectra  as  a function  of  g-value. 
All  spectra  measured  at  30  K,  and  240  GHz.  Each  spectrum  was  signal 
averaged  (10-20  times),  which  takes  3 to  7 hours. 
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Figure  4-9  Effective  hyperfme  splittings  obtained  from  the  fitting  of  the  experimental 
ENDOR  spectra  shown  in  figure  4-8. 
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Bacteriopheophytin  a Radical  Anion 

BPheo  a radical  anion  was  generated  chemically,  as  described  in  Chapter  3.  This 
radical  was  measured  with  EPR,  ENDOR  and  ESEEM  at  X-band,  as  well  as  EPR  and 
ENDOR  at  W-band  and  at  240  GHz.  The  cw-EPR  at  X-Band  shows  partially  resolved 
hyperfine  structure  at  room  temperature.  Figure  4-10  shows  the  experimental  and  the 


Magnetic  Field  (G) 

Figure  4-10  Cw-EPR  of  BPheo  a radical  anion  at  298  K and  9.42581  GHz  (X-band),  with 
1.0  mW  microwave  power,  10ms  time  constant,  20  ms  conversion  time,  and 
0.7  G field  modulation  at  100  kHz.  The  experimental  data  (solid)  were 
simulated  (dotted)  to  determine  three  unique  isotropic  hyperfine  coupling 
constants:  3 x Amethine  = 2.5  G,  6 x Amethyi  = 3.0  G,  and  2 x AN  = 2.3  G. 

simulated  X-band  spectra.  This  experiment  was  done  at  298  K,  1.0  mW  microwave 

power,  0.5  G field  modulation  at  100  kHz,  10  ms  time  constant,  20  ms  conversion  time, 

and  a microwave  frequency  of  9.42581  GHz.  The  experimental  data  (solid)  was  fitted 

using  the  DOUBLET  software  to  determine  giS0  = 2.003263  ± 6x1  O'6.  The  simulation 

(dotted)  used  averaged  contributions  from  the  methine  bridge  protons  (a,  (3,  8),  methyl 

protons  (la  and  5a),  and  nitrogens  22  and  24  (see  figure  1-3).  The  hyperfine  coupling 
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constants  used  for  the  fitting  are  as  follows:  Amethine  = 2.5  G (7.0  MHz),  Amethyi  = 3.0  G 
(8.4  MHz),  and  AN  = 2.3  G (6.4  MHz).  The  linewidth  used  in  the  simulation  was  0.92  G. 
These  hyperfine  couplings  are  in  agreement  with  previously  published  results,  which  are 
collected  with  the  experimental  results  in  table  4-1  [56], 

Table  4-1  Isotropic  hyperfine  coupling  constants  for  BChl  a ' and  BPheo  a ~ determined 
with  X-band  EPR. 


Bacterial 

Coupling 

# atoms 

hfc  from  EPR 

hfc  from  ENDOR 

hfc  from 

Pigment 

Assignment 

(Gauss)  [56] 

(Gauss)  [56] 

experiment 

(Gauss) 

BChl  a ’ 

Amethine 

2 

2.5 

2.5 

2.4 

Amethine 

1 

2.0 

2.0 

2.4 

an 

2 

2.2 

2.2 

2.2 

Amethyi 

6 

3.2 

3.1,  3.3 

3.0 

BPheo  a ' 

Amethine 

3 

2.5 

2.5 

2.5 

An 

2 

2.3 

2.3 

2.3 

Amethyi 

6 

3.0 

3.0,  3.0 

3.0 

Aunre  solved 

4 

<1 

0.1,  0.6 

0.92  G 
(linewidth) 

In  addition  to  the  cw  experiments  done  at  X-band,  pulsed  experiments  were  done 
in  an  attempt  to  estimate  the  hyperfine  splittings  of  the  coupled  nuclei.  Figure  4-11 
shows  the  results  of  a 3-pulse  ESEEM  experiment  done  at  23  K.  The  pulse  separation 
between  the  second  and  third  pulses  was  increased  in  8 ns  steps  from  100  ns  to  16,376  ns. 
This  ESEEM  spectrum  was  baseline  corrected  for  the  Ti  decay  and  Fourier  transformed 
using  the  Bruker  software  on  the  E580  spectrometer.  The  Fourier  transform,  shown  in 
figure  4-12,  shows  a somewhat  weak  proton  coupling  at  14.7  MHz,  as  well  as  peaks  at 
1.77  MHz  and  2.35  MHz,  which  result  from  the  single  quantum  transitions  between 
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quadrupole  levels  of  the  nitrogen  nuclei  (1=1).  Isotropic  hyperfine  coupling  constants 
may  be  estimated  from  the  frequencies  of  these  nitrogen  transitions,  but  this  data  is 
known  and  was  not  extracted  here. 


Figure  4-11  X-band  3-pulse  ESEEM  on  BPheo  a radical  anion  at  23K.  The  T i decay 
curve  was  fitted  and  subtracted  from  the  ESEEM  spectrum  to  generate  this 
plot.  The  microwave  frequency  was  9.6854  GHz,  with  the  magnetic  field  at 
3452  G,  with  16  ns  pulses,  x = 120  ns,  and  T incremented  in  4 ns  steps. 

Pulsed  ENDOR  experiments  were  done  at  W-Band  (95  GHz)  on  a Bruker  e680 
spectrometer  at  the  NHMFL  in  Tallahassee,  FL.  These  experiments  showed  slight  g- 
tensor  resolution  enhancement  over  X-band,  as  seen  in  figure  4-13.  This  spectrum  is  a 
plot  of  the  Hahn  echo  intensity  as  a function  of  the  external  field.  The  small  shoulder  on 
the  low  field  side  of  the  signal  belongs  to  g±  = 2.007766  while  the  rest  of  the  signal 
belongs  to  g(  = 2.005667.  These  g-values  are  likely  inaccurate  due  to  the  uncertainty  in 
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Figure  4-12  Fast  Fourier  transform  of  the  time  domain  3-pulse  ESEEM  data  in  figure  4- 
11.  The  spectrum  shown  in  figure  4-11  was  zero-filled  to  2046  points,  then 
smoothed  with  an  RC  filter  before  FFT  to  generate  this  plot. 

the  static  magnetic  field,  but  the  anisotropy  is  consistent  with  the  expected  high  field  g- 

values,  since  Ag  = 2.09x10°.  ENDOR  experiments  were  done  at  the  fields  indicated  in 

figure  4-13.  Figure  4-14  shows  the  Davies  ENDOR,  which  employed  a long,  320  ns 

saturating  microwave  pulse,  a 12  microsecond  RF  pulse,  a short  180  ns  microwave  pulse, 

and  a 236  ns  refocusing  pulse.  The  horizontal  scale  on  the  figure  is  the  difference 

between  the  RF  frequency  of  the  pulse  and  the  proton’s  Larmor  frequency.  The  central 

feature  of  the  spectra,  at  0 MHz  offset  frequency,  results  from  the  matrix  protons  in  the 

sample  that  show  no  hyperfine  splitting.  A small  splitting  of  about  2 MHz,  is  apparent, 

yet  unresolved,  and  appears  as  shoulders  on  the  matrix  peak.  ENDOR  transitions  are  also 
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Figure  4-13  W-band  echo  detected  field  sweep  of  BPheo  a radical  anion,  measured  at 
15K  and  94.08275  GHz.  The  microwave  pulses  were  72  ns  and  120  ns  in 
length,  separated  by  300  ns.  The  Ag  between  33508  G (g  = 2.00777)  and 
33543  G (g  = 2.00567)  is  2.09xl0'3,  which  fits  the  estimated  g-tensor 
anisotropy  obtained  at  high  field  EPR  measurements. 

visible  around  + 4.0  MHz,  at  A = 8.0  MHz.  Despite  the  presence  of  this  ENDOR 

intensity  around  + 4 MHz,  the  lack  of  resolution  prevents  the  accurate  determination  of 

the  hyperfine  coupling.  Also,  several  hyperfine  couplings  are  expected  to  be  observed  in 

this  region,  so  this  intensity  is  most  likely  a superposition  of  several  ENDOR  lines.  If 

orientation  selection  has  sufficiently  resolved  the  two  apparent  g-components,  the 

convolution  of  lines  has  obscured  the  anisotropic  information.  These  results  indicate  the 

need  for  higher  fields  to  resolve  the  anisotropic  hyperfine  couplings. 
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Figure  4-14  Davies  ENDOR  of  BPheo  a radical  anion  taken  at  15K  and  94.077  GHz  + 
0.005  GHz.  This  experiment  employed  a 320  ns  saturating  microwave  pulse, 
a 12  p,s  RF  pulse,  a 180  ns  microwave  pulse,  and  a 236  ns  refocusing  pulse. 
Each  spectrum  is  a result  of  20-40  averages,  7-10  hours/spectrum. 

BPheo  a radical  anion  was  measured  at  high-field  using  the  transmission  cw-EPR 
setup  described  in  chapter  3.  Figure  4-15  shows  the  high  field  EPR  spectrum  (solid)  as 
well  as  a simulated  powder  pattern  (dotted).  The  experiment  was  done  at  10  K with 
~1  mW  microwave  power,  4 Gauss  field  modulation  at  1 kHz,  and  a microwave 
frequency  of  326.577  GHz.  The  sample  was  accompanied  by  a phosphorus  doped  silicon 
field  calibration  standard.  The  doublet  lines  at  11.6750  T result  from  the  phosphorus 
signal  (1=  Vi)  in  the  standard,  while  the  feature  at  1 1.6400  T is  of  unknown  origin.  The 
simulation  of  the  EPR  data  was  performed  using  DOUBLET  software  that  calculates 
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Figure  4-15  High  field  EPR  of  BPheo  a radical  anion  in  THF  at  10  K and  326.577  GHz, 
with  lmW  microwave  power,  and  6 G field  modulation  at  ~1  kHz,.  The 
experiment  data  (solid)  was  simulated  (dotted),  employing  anisotropic 
linewidths  for  the  components  of  g:  Ax  = 7 G,  Ay  = 7 G,  Az  = 15  G. 

field  values  based  on  the  work  of  Michael  Zemer  [80].  The  phosphorus  spectral 

contribution  was  fitted  and  the  g-value  obtained  from  that  fitting,  g(Pexp),  was  used  to 

correct  the  g-values  for  the  powder  pattern.  The  powder  pattern  was  fitted  and  the 

experimental  g-values  from  that  fitting  were  adjusted  by  the  correction  factor,  which  is 

the  ratio  of  the  accepted  phosphorus  g-value,  g(P)  = 1.998525  + 5xl0'5,  and  the 

experimental  g(Pexp).  The  fitting  and  the  field  correction  indicate  the  following  g-values: 

gx  = 2.004459  ± 5x1  O'5,  gy  = 2.003249  ± 5xl0'5,  and  gz  = 2.002299  ± 5xlO'5.  The 

average  of  these  principal  components  of  g give  the  isotropic  constant,  gave  = 2.00334  + 

5xl0‘5.  The  estimated  error  in  these  values  is  derived  from  the  minimum  amount  of 

variability  in  the  fitted  g-values  while  maintaining  a good  fit  to  the  experimental 
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spectrum.  The  discrepancies  between  the  calculated  and  experimental  intensities  may 
result  from  incorrect  polarization  of  the  microwaves  throughout  the  sample.  Without  a 
cavity,  the  polarization  of  the  microwave  radiation  may  not  be  consistent  throughout  the 
sample. 

This  sample  was  also  measured  on  the  superheterodyne  setup  at  NHMFL,  which 
has  ENDOR  capabilities.  These  cw-ENDOR  experiments  detect  changes  in  the  EPR 
intensity  as  a function  of  the  RF  frequency.  The  EPR  transition  is  saturated,  shown  in 
figure  4-16,  then  ENDOR  experiments  were  done  at  12  different  field  positions  across  the 
powder  spectrum.  These  ENDOR  spectra  are  shown  in  figures  4-17  to  4-19  and  show 
hyperfine  interactions  that  are  field  dependent.  These  ENDOR  spectra  are  collected  in  a 
stack  plot  in  figure  4-20,  where  the  ENDOR  frequency  seems  to  be  dependent  on  the  g- 
value. 


Figure  4-16  Saturated  EPR  absorption  of  BPheo  a radical  anion  at  40  K and  239.57  GHz, 
with  ~1  mW  microwave  power  and  6 G field  modulation  at  1 kHz.  The 
underlying  backgound  signal  (below  -13  mT  and  above  5 mT)  is  due  to 
manganese  in  the  sample  holder.  Rapid  passage  and  saturation  affects  lead  to 
the  absorptive  appearance  of  this  powder  pattern. 
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Figure  4-17  High  field  cw-ENDOR  on  BPheo  a anion,  measured  on  gx  at  g = 2.00415 
and  10  K.  The  experimental  spectrum  (solid)  was  simulated  using  the 
GENDOR  software  (dotted)  to  estimate  the  effective  hyperfine  splittings, 
shown  in  figure  4-21,  and  the  anisotropic  hyperfine  tensors  discussed  in 
chapter  5. 


Figure  4-18  High  field  cw-ENDOR  of  BPheo  a anion,  measured  on  gy  at  g = 2.00338  and 
10  K.  The  experimental  spectrum  (solid)  was  simulated  using  the  GENDOR 
software  (dotted)  to  estimate  the  effective  hyperfine  splittings,  shown  in  figure 
4-21,  and  the  anisotropic  hyperfine  tensors  discussed  in  chapter  5. 
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Figure  4-19  High  field  cw-ENDOR  of  BPheo  a anion,  measured  on  gz  at  g = 2.00228  and 
10  K.  The  experimental  spectrum  (solid)  was  simulated  using  the  GENDOR 
software  (dotted)  to  estimate  the  effective  hyperfine  splittings,  shown  in  figure 
4-21,  and  the  anisotropic  hyperfine  tensors  discussed  in  chapter  5. 
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Figure  4-20  Stack  plot  of  cw-ENDOR  spectra  on  BPheo  a radical  anion  at  10  K and 

239.57  GHz.  Each  spectrum  was  signal  averaged  10-30  times,  which  requires 
4-9  hours/spectrum. 

The  position  of  the  ENDOR  transitions  leads  to  effective  hyperfine  splittings. 


which  are  presented  in  figure  4-21. 
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Figure  4-21  Effective  hyperfme  splittings  for  BPheo  a radical  anion  obtained  through 

simulation  of  the  orientationally  selective  ENDOR  spectra  shown  in  figure  4- 

20. 
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Bacteriochlorophyll  a Radical  Cation 

The  BChl  a radical  cation  was  investigated  using  pulsed  W-band  EPR  and  ENDOR 
techniques.  The  EPR  signal  was  initially  detected  using  X-band,  which  shows  little 
hyperfine  resolution,  shown  in  figure  4-22.  This  spectrum  is  shown  with  a DOUBLET 
simulation,  which  determined  giso  = 2.003308  + 5xl0'5.  The  simulation  includes  isotropic 
hyperfine  interactions  that  were  determined  through  ENDOR  and  TRIPLE  resonance 
studies  [67,109],  Higher  fields  were  employed  in  attempts  to  resolve  the  anisotropic 
interactions.  The  W-band  EPR  signal  was  detected  using  a Hahn  echo  pulse  sequence. 
This  experiment  uses  a 72  ns  pulse,  followed  by  a 120  ns  pulse.  Figure  4-23  plots  the 


Magnetic  Field  (Gauss) 

Figure  4-22  X-band  cw-EPR  spectrum  of  BChl  a radical  cation  at  298  K and  9.440197 
GHz  with  2 mW  microwave  power,  10  ms  time  constant,  20  ms  conversion 
time,  and  1 G field  modulation  at  100  kHz.  The  experimental  data  (solid)  was 
fitted  using  DOUBLET,  which  employed  a Lorentzian  lineshape  with  3.9  G 
linewidth  and  the  following  hyperfine  couplings:  4 x Abeta=  4 G,  3 x Amethyi-i 
= 1.76  G,  and  3 x Amethyi-5  = 3.4  G. 
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intensity  of  the  spin  echo  as  a function  of  magnetic  field  at  W-band.  This  EPR  line 
appears  to  be  a single,  featureless  line  with  little  g-resolution.  This  result  is  not 
surprising  considering  the  g-tensor  anisotropy  for  this  cation  is  known  to  be  less  than  the 
corresponding  anion.  The  spin-spin  (T2)  and  spin  lattice  (TO  relaxation  times  were 
determined  to  be  475  ns  and  6 p.s,  respectively.  T2  was  measured  by  generating  a spin 
echo,  then  increasing  the  time  between  the  second  microwave  pulse  and  the  detection,  x. 
This  separation,  t,  is  plotted  against  the  echo  intensity.  This  experiment  uses  two  120  ns 
microwave  pulses,  initially  separated  by  400  ns,  which  is  increased  in  8 ns  steps.  This 
experiment  detects  the  diffusion  of  the  spin  packets,  which  is  seen  as  a 


Magnetic  Field  (Gauss) 

Figure  4-23  Echo  detected  field  sweep  on  BChl  a radical  cation  measured  at  15  K and 

94.2862  GElz.  Two  80  ns  pulses  were  separated  by  240  ns,  followed  by  a 180 
ns  pulse  at  604  ns.  The  arrows  indicate  the  field  values  for  the  Davies 
ENDOR  spectra. 
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reduction  in  the  echo  intensity  as  x is  increased.  T2  was  found  to  be  475  ns.  The  Ti 
measurement  was  estimated  using  a saturation  recovery  experiment,  which  polarizes  the 
spins  by  180°  microwave  pulse,  followed  by  a long  delay  time  (on  the  order  of  Tj),  which 
is  increased  in  32  ns  steps.  The  magnetization  is  then  projected  onto  the  xy-plane  by  a 
Hahn  echo  sequence.  Thus,  the  magnitude  of  the  longitudinal  magnetization  vector  will 
determine  the  intensity  of  the  echo.  The  echo  intensity  is  monitored  as  the  delay  between 
the  long  saturating  pulse  and  the  Hahn  sequence  are  increased.  The  echo  intensity  grows 
with  a time  constant,  Ti  = 6 ps,  for  this  sample. 

Mims  and  Davies  ENDOR  experiments  were  performed  at  the  field  positions 
indicated  in  figure  4-23.  These  ENDOR  spectra  are  shown  in  figures  4-24  and  4-25.  The 
Davies  ENDOR  experiments  use  a 184  ns  microwave  pulse,  followed  by  an  18  ps  RF 


Figure  4-24  Davies  ENDOR  on  BChl  a radical  cation  taken  at  15K  and  94.1415  GHz, 

which  employed  a 180  ns  n pulse,  1 1 ps  RF  pulse,  and  two  160  ns  microwave 
pulses.  The  shot  repetition  time  was  8 ms. 
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pulse,  followed  by  a Hahn  echo  pulse  sequence  (92ns  and  160  ns).  Figure  4-24  shows  the 
Davies  spectra  taken  at  different  field  positions  on  the  EPR  envelope.  The  hyperfine 
couplings  show  little  variation  between  these  spectra,  which  indicates  that  orientational 
selection  may  not  have  been  achieved.  However,  the  spectra  taken  at  the  low  and  high 
field  edges  of  the  EPR  envelope  show  an  additional  hyperfine  coupling  which  is  not 
visible  in  the  33585  G spectrum.  The  experimental  values  for  these  couplings  are 
tabulated  in  table  4-2.  The  Mims  ENDOR  spectra  show  similar  spectra  at  the  different 
field  positions,  however,  the  differences  are  subtle,  and  unfortunately,  not  resolved 
enough  to  make  definitive  assignments. 


v - vH  (MHz) 

Figure  4-25  Mims  ENDOR  on  BChl  a radical  cation  taken  at  15K  and  94.2862  GHz  with 
two  80  ns  microwave  pulses  separated  by  240  ns,  a 12  ps  RF  pulse,  then  a 
third  80  ns  microwave  pulse.  The  shot  repetition  time  was  8 ms. 
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Table  4-2  Davies  and  Mims  ENDOR  data  for  BChl  a radical  cation 


Field  (G) 
Davies 

Field  (G) 
Mims 

hfc  # 

A +0.25  MHz  A + 0.25  MHz 

Davies  Couplings  Mims  Couplings 

33530 

33570 

1 

1.18 

1.16 

2 

4.93 

5.12 

3 

9.39 

- 

33540 

33580 

1 

1.49 

1.12 

2 

5.17 

5.1 

3 

9.86 

- 

33585 

33604 

1 

- 

1.11 

2 

4.93 

5.09 

3 

9.78 

- 

33615 

33620 

1 

1.48 

1.17 

2 

5.09 

4.63 

3 

9.24 

- 

33625 

1 

- 

1.22 

2 

- 

4.69 

Bacteriopheophytin  a Radical  Cation 

The  BPheo  a radical  cation  was  successfully  studied  with  X-band  and  W-band 
EPR  and  ENDOR,  while  complications  limited  the  success  of  the  high  field  work.  The 
X-band  EPR  results  in  a featureless  line  at  room  temperature  and  under  cryogenic 
conditions,  with  a linewidth  of  9 G (not  shown).  W-band  experiments  showed  little 
enhancement  in  g-resolution,  as  seen  in  figure  4-26  by  the  lack  of  resolution  in  the 
principal  g-components.  The  bulk  of  this  signal  results  from  the  BPheo  a cation,  while 
the  sharp  feature  at  33668  G could  be  another  radical  species  formed  during  the 
oxidation.  The  33668  G feature  is  most  likely  due  to  a different  radical  with  a much 
longer  T2,  since  it  has  a much  sharper  linewidth  than  the  BPheo  a radical  signal.  This 
sharp  radical  signal  has  g = 2.002964,  while  the  low  field  signal  is  closer  to  2.005943  + 
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0.003.  Unfortunately,  the  lack  of  a field  calibration  standard  prevents  accurate  g-value 
measurements  at  this  frequency.  ENDOR  spectra  were  taken  at  the  indicated  points  on 
the  spin  echo  envelope.  The  Davies  spectra,  shown  in  figure  4-27,  show  little  anisotropy 
with  respect  to  the  magnetic  field.  The  measured  hyperfine  coupling  values  from  the  W- 
band  experiments  are  listed  in  table  4-5.  The  Mims  ENDOR  spectra  show  slight 
differences  in  the  range  of  A = 6 MHz.  The  high  field  spectrum  in  figure  4-28  shows 
broad  wings  with  A = 14-15  MHz,  which  are  absent  in  the  low  field  spectrum.  Due  to 
the  apparent  lack  of  resolution  in  g,  these  experiments  were  not  pursued.  High  field  EPR 
was  attempted  but  line-broadening  affects  were  observed  above  W-band,  which 
prevented  the  measurement  of  the  principal  components  of  g. 


Magnetic  Field  (Gauss) 

Figure  4-26  Echo  detected  field  sweep  of  BPheo  a radical  cation  at  20  K and  94.3859 

GHz.  This  experiment  employed  two  microwave  pulses  of  length,  192  ns  and 
260  ns  separated  by  800  ns.  The  shot  repetition  time  was  6 ms. 
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Figure  4-27  Davies  ENDOR  of  BPheo  a radical  cation  at  20  K and  94.0827  GHz,  which 
employed  a 200  ns  n pulse,  12  ps  RF  pulse,  and  two  190  ns  microwave 
pulses.  The  shot  repetition  time  was  8 ms. 


Figure  4-28  Mims  ENDOR  of  BPheo  a radical  cation  at  20  K and  94.0730  GHz,  which 
employed  three  80  ns  microwave  pulses  with  the  second  and  third  pulses 
separated  by  a 12  ps  RF  pulse.  The  shot  repetition  time  was  10  ms. 
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Table  4-3  BPheo  a radical  cation  pulsed  W-band  ENDOR  data 


ENDOR 

Experiment 

Field 

Ai 

A2 

a3 

A4 

Davies 

33618 

1.6  MHz 

3.3  MHz 

5.2  MHz 

8.7  MHz 

33668 

1.6  MHz 

3.1  MHz 

5.2  MHz 

7.7  MHz 

Mims 

33618 

1.3  MHz 

5.2  MHz 

33668 

1.3  MHz 

5.2  MHz 

Chlorophyll  a Radical  Anion 

The  Chi  a radical  anion  was  generated  using  potassium  metal  in  tetrahydrofuran 
(THF)  using  high  vacuum  techniques  described  in  chapter  3.  This  radical  was  examined 


with  EPR  at  9.5  GHz  and  330  GHz.  X-Band  EPR  spectra  consist  of  a single 
inhomogeneously  broadened  line,  with  an  isotropic  g-value  of  2.002667  + 6xl0'6.  Figure 
4-29  shows  the  X-band  spectrum  taken  at  298  K,  with  a 100  kHz  field  modulation  of  2 G, 


20.48  ms  time  constant,  40.96  ms  conversion  time,  and  2 mW  microwave  power.  This 
radical  signal  has  a linewidth  of  8 G at  room  temperature.  The  experimental  data  was 
simulated  using  DOUBLET  to  obtain  the  g-value,  as  well  as  the  indication  of  unresolved 


hyperfine  interactions.  The  fit  (dotted  line)  in  figure  4-29  uses  hyperfine  contributions 
from  three  different  types  of  protons  (A  = 1.1,  2.7,  2.9  G).  These  contributions  are 


estimates  since  the  splitting  is  unresolved  under  these  conditions.  The  lineshape  in  this 
simulation  is  assumed  to  be  90  percent  Lorentzian.  EPR  of  this  radical  was  also 


measured  at  higher  fields  in  frozen  solution,  shown  in  figure  4-30.  This  experiment  was 
performed  on  the  transient  EPR  setup  at  the  NHMFL,  described  in  chapter  3.  The  g- 
tensor  achieved  full  resolution  above  280  GHz,  with  values:  gx  = 2.004074  ± 5xl0"5,  gy  = 
2.002884  + 5x10 5,  gz=  2.002024  + 8xl0‘5,  which  lead  to  an  average  value  of  gave  = 
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Figure  4-29  EPR  spectrum  of  Chi  a radical  anion  at  298  K and  9.423040  GHz  with  2 
mW  microwave  power  and  2 G field  modulation  at  100  kHz.  The 
experimental  data  (solid)  was  simulated  (dotted)  assuming  a Lorentzian 
lineshape  at  gjso=  2.002667  + 6x10  3 with  three  different  proton  hyperfine 
interactions:  A|  = 1.1G,A2  = 2.7,  and  A3  = 2.9  G. 

2.002994  + 8xl0"5.  The  g-values  were  obtained  through  fitting  the  data  with  the 

DOUBLET  software.  The  magnetic  field  was  calibrated  using  a phosphorus  doped 

silicon  standard  with  known  g-value.  The  field  correction  procedure  was  described  in  the 

bacteriopheophytin  a radical  anion  high  field  EPR  results.  This  spectrum  was  fitted  using 

anisotropic  linewidth  approximations  where,  Ax  = 10  G,  Ay  = 7 G,  and  Az  = 15  G.  The 

poor  fit  on  gz  introduces  some  error  in  this  data  point,  and  the  sources  of  this  poor  fit  will 


be  addressed  in  the  discussion. 


74 


Figure  4-30  EPR  of  Chi  a radical  anion  with  phosphorus  standard  at  15  K and  326.13 
GHz  with  4 G field  modulation,  and  ~0.9  mW  microwave  power.  The 
experimental  spectrum  (solid)  was  simulated  (dotted)  using  anisotropic 
linewidths  given  by:  Ax  = 10  G,  Ay  = 7 G,  Az  = 15  G.. 

Chlorophyll  b 

Chi  b radical  anion  was  studied  with  EPR  at  X-band  and  at  330  GHz.  The  X- 
band  spectrum  at  298  K shows  a single  EPR  line  with  giso=  2.002846  + 6x10  6.  Figure  4- 
31  shows  the  X-band  data  with  a simulation  that  accounts  for  three  predicted  hyperfine 
couplings.  The  parameters  for  this  experiment  are:  2 G field  modulation  at  100  kHz, 
20.48  ms  time  constant,  40.96  ms  conversion  time,  and  2 mW  microwave  power. 
Unresolved  hyperfine  couplings  are  present  and  are  possibly  broadened  by  concentration 
effects.  However,  these  hyperfine  couplings  must  be  on  the  order  of  the  linewidth,  so 
they  should  hardly  be  resolved  under  these  conditions.  The  solid  line  in  figure  4-31  is 
experimental  data  and  the  dotted  line  is  a simulation  obtained  with  the  DOUBLET 
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Figure  4-31  X-band  EPR  of  Chi  b radical  anion  at  298K  and  9.426530  GHz,  with  1 mW 
microwave  power,  10  ms  time  constant,  20  ms  conversion  time,  and  1 G field 
modulation  at  100  kHz.  The  experimental  data  (solid)  was  simulated  (dotted) 
using  the  DOUBLET  software  with  the  following  hyperfine  coupling 
constants:  3 x Ai  = 1.0  G,  3 x A2  = 2.7  G,  and  3 x A3  = 3.1  G. 

software.  The  simulation  approximates  a Lorentzian  shape,  with  isotropic  hydrogen  and 

nitrogen  couplings,  A < 3 G. 

The  Chi  b radical  anion  was  measured  using  high  field,  cw-EPR  on  the 
transmission  instrument  described  in  chapter  3.  This  experiment  applies  326.13  GHz 
microwave  radiation  at  magnetic  fields  around  1 1 T.  These  experimental  conditions  are 
required  to  fully  resolve  the  g-tensor  into  the  three  principle  components.  The  high  field 
spectrum  is  shown  in  figure  4-32,  with  gx=  2.003964  + 5x10  5,  gy  = 2.002779  ± 5xlO"5, 
and  gz  = 2.001896  + 7xl0';\  which  lead  to  an  average  value  of  gave  = 2.002880  + 7x10  5. 
The  radical  sample  is  accompanied  by  a phosphorus  doped  silicon,  field  calibration 
standard.  The  standard  is  the  source  of  the  weak  doublet  up  field  from  the  powder 
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pattern.  The  field  correction  procedure  has  been  discussed  in  the  BPheo  a radical  anion 
results.  The  experimental  g-values  are  shown  in  table  4-4. 


Magnetic  Field  (Tesla) 

Figure  4-32  High  field  EPR  spectrum  of  Chi  b radical  anion  at  15K  and  326.13  GHz, 
with  5 G field  modulation  and  ~0.9  mW  microwave  power.  The  two  weak, 
high  field  lines  are  the  phosphorus  field  calibration  standard.  The 
experimental  spectrum  (solid)  was  simulated  (dotted)  using  anisotropic 
linewidths  given  by:  Ax  = 10.5  G,  Ay  = 7.3  G,  Az  = 17  G.. 
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Table  4-4  Radical  anion  g-tensor  data.  These  measurements  were  taken  on  the 

transmission  setup  and  correspond  to  different  samples,  which  were  prepared 
identically,  according  to  the  procedures  in  chapter  3. 


Sample 

Sxx 

gyy 

§ZZ 

(§xx'8yy) 

(8yy”Szz) 

Microwave 

Frequency 

(GHz) 

Bacteriochlorophyll  a ' 

2.004257 a 

2.003150“ 

2.002513“ 

0.001107 

0.000637 

326.81 

2.004257 a 

2.003123“ 

2.002333 b 

0.001134 

0.000790 

327.98 

2.004237 a 

2.003112“ 

2.00205  lc 

0.001125 

0.001061 

326.81 

2.004567 

2.003207 

2.002256 

0.001360 

0.00951 

DFT 

Bacteriopheophytin  a ' 

2.004459 a 

2.003249“ 

2.002299“ 

0.001210 

0.000950 

326.13 

2.004423 a 

2.003 144 c 

2.002274 f 

0.001279 

0.000870 

326.13 

2.004504 b 

2.003294 b 

2.002344 e 

0.001210 

0.000950 

326.58 

2.003991 

2.002945 

2.002176 

0.001047 

0.000768 

DFT 

Chlorophyll  a ' 

2.004012“ 

2.002912“ 

2.002027 d 

0.00110 

0.000885 

326.11 

2.004074a 

2.002884“ 

2.002024 d 

0.00119 

0.000860 

326.36 

2.004079“ 

2.002869“ 

2.002239 f 

0.00121 

0.00063 

326.16 

Chlorophyll  b ' 

2.003964 b 

2.002779 b 

2.001 896 c 

0.001185 

0.000883 

326.13 

The  superscripts  a-f  denote  the  experimental  error  in  that  measurement:  a = 5xlO'5,  b = 6xlO'5, 
c = 7x10 5,  d = 8xl0‘5,  e = 9xl0'5,  f = 1x10 4.  These  values  were  derived  from  the  quality  of  the  fitting 
to  the  hf-EPR  data. 

CHAPTER  5 
DISCUSSION 

Bacteriochlorophyll  a Radical  Anion 

At  210  GHz  the  anisotropy  of  the  Zeeman  splitting  exceeds  the  hyperfine 
linewidth,  allowing  accurate  g-tensor  determination.  The  high-field  EPR  spectrum 
shown  in  figure  4-3  shows  a fully  resolved  powder  pattern  with  rhombic  symmetry.  The 
following  discussion  will  require  the  use  of  a reference  frame  for  the  g-tensor  orientation. 
The  molecular  axes  will  be  defined  X,  Y,  and  Z,  where  X forms  a vector  between  Nn  and 
Niv,  Y forms  a vector  between  Ni  and  Nm,  and  Z is  normal  to  the  XY  plane  containing 
all  four  nitrogens.  The  gz=  2.002513  + 5x1  O'5  component  closely  resembles  the  free 
electron  g-value,  ge=  2.0023193.  This  value  is  in  good  agreement  with  density 
functional  theory  (DFT)  predictions  using  a double  zeta  basis,  which  gives  gz  = 
2.0022558,  oriented  84°  to  the  XY  plane.  The  IIRC  calculation  used  an  unrestricted 
Hartree-Fock  double  zeta  basis,  which  describes  each  atomic  orbital  with  two  functionals 
(dz)  that  are  functions  of  the  electron  density  and  the  gradient  of  this  density  over  the 
entire  molecule.  The  gradient  correction  to  the  exchange  functional  is  “revised  PBE”  and 
the  correlation  functional  is  “PBE.”  This  predicted  orientation  of  gz  is  in  line  with 
generally  accepted  theory  regarding  the  nature  of  the  electron  density  on  a planar,  organic 
7t-radical.  According  to  Stone,  the  gz  axis  for  doublet  ^-radicals  should  lie  perpendicular 
to  the  plane  of  conjugation  [76].  In  this  normal  orientation,  the  pz  orbital  cannot  mix  with 
heavy  atoms  or  orbitals  in  the  plane  of  the  conjugation  to  obtain  orbital  angular 
momentum.  This  orientation  allows  it  to  maintain  a value  close  to  that  of  the  free 
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electron,  as  shown  in  equation  2-3.  Relativistic  mass  effects  and  gauge  corrections  are 
proposed  to  explain  the  observed  behavior  of  gz<  ge  [41,75,79,80],  However,  the 
deviation  of  the  measured  gz  value  from  ge  could  also  result  from  an  angle  slightly  less 
than  90°  to  the  conjugation  plane,  as  was  determined  by  Bratt  et  al.  for  the  Chi  a radical 
cation  [74], 

In  order  to  fit  the  powder  pattern,  anisotropic  linewidths  for  the  principal 
components  of  g were  employed.  These  linewidth  affects  indicate  a broadening 
phenomenon  such  as  anisotropic  hyperfine  interaction  or  environmental  interactions  that 
selectively  broaden  gz.  Burghaus  et  al.  have  observed  selective  broadening  effects  in 
quinone  radicals  where  hydrogen-bonding  interactions  shift  the  components  of  g oriented 
toward  the  interacting  solvent  molecule  up  to  g + lxlO"4  [113],  Their  study  observed 
dipolar  broadening  of  certain  lines  of  a powder  pattern  depending  on  the  hydrogen 
bonding  ability  of  the  solvent  with  the  quinone  radical.  A similar  phenomenon  may  be 
responsible  for  the  broadened  gz  component  of  the  BChl  a radical  anion  EPR  spectrum. 
The  THF  solvent  molecules  contain  an  oxygen  atom  which  can  serve  as  a fifth  ligand  to 
the  central  Mg  atom  in  BChl  a.  Although  hydrogen  bonding  would  be  expected  to  have  a 
stronger  effect  on  the  MO  of  the  radical,  the  potential  for  electrostatic  interaction  exists. 
This  prediction  is  justified  by  the  localization  of  electron  density  in  the  Z dimension  of 
the  radical,  which  is  the  very  same  dimension  that  has  the  potential  for  intermolecular 
interactions  with  Mg.  It  should  be  mentioned,  that  the  concentration  of  spin  density  on 
the  oxygen  atoms  of  quinone  radicals  is  much  higher  than  on  the  magnesium  atom  of 
BChl  a radical  anion.  This  behavior  should  weaken  the  effect  of  intermolecular 


interactions  on  the  g-tensor. 
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The  remaining  components  of  the  g-tensor  were  found  to  be  gx  = 2.004257 
+ 5x10 3,  gy  = 2.003150  + 5xl0'5,  which  deviate  from  gz.  The  spin-orbit 
interactions  that  lead  to  this  deviation  mainly  result  from  heavy  atoms,  such  as  nitrogen 
and  oxygen.  The  nitrogen  atoms  of  rings  II  and  IV  are  directly  linked  to  the  conjugated 
system,  which  should  lead  to  a higher  concentration  of  spin  density  than  Ni  and  Nm.  This 
assumption  would  indicate  a larger  degree  of  spin-orbit  interaction  along  the  X-axis  than 
on  the  Y-axis  of  the  radical.  Only  two  oxygen  atoms,  the  ring  I acetyl  group  and  the  ring 
V carbonyl  group,  are  expected  to  interact  with  the  spin  density  of  this  radical.  It  has 
been  predicted  that  the  acetyl  group’s  orientation  can  affect  the  values  of  the  g-tensor  and 
it's  orientation  for  BChl  a+  [33,73]  and  Chi  a [74],  Plato  et  al.  predict  that  by  rotation  of 
the  acetyl  group,  relative  to  the  conjugation  plane,  gz  and  gy  become  nearly  degenerate  at 
an  angle  of  + 45°  [33].  As  this  group  is  rotated  out  of  the  molecular  plane,  the  pz  orbitals 


Figure  5-1  Proposed  orientation  of  g on  the  BChl  a radical  anion.  The  smallest  spin-orbit 
interactions  lead  to  gz  oriented  perpendicular  to  the  plane  of  the  molecule, 
while  the  highest  degree  of  spin-orbit  interactions  are  found  on  rings  I and  III. 
The  nitrogen  atoms  on  rings  II  and  IV  also  interact  with  the  spin  density, 
bringing  the  value  of  gy  away  from  gz. 
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can  mix  with  the  lone  pairs  on  the  oxygen.  This  interaction  increases  gz  and  decreases  gy, 
which  results  in  near  axial  symmetry.  The  g-tensor  will  be  influenced  by  the  amount  of 
spin  density  on  the  acetyl  group,  which  may  differ  between  the  cation  and  the  anion. 

Since  powder  EPR  experiments  provide  limited  information  on  the  orientation  of  the  g- 
tensor,  theoretical  models  were  employed. 

Using  the  revPB  BLYP  level  of  theory  and  employing  a double-zeta  basis  set, 
DFT  predicts  the  following  g-values:  gx=  2.004567  and  gy  = 2.003207,  which  differ  from 
experiment  by  2xl0"4  for  gx  and  -2.3xl0"5  for  gy.  The  calculation  of  gx  predicted  an 
orientation  of  89°  from  X and  105°  from  Y,  while  gy  was  predicted  to  be  19°  from  X and 
74°  from  Y.  This  orientation  of  the  g-tensor  on  the  molecular  axes  is  depicted  in  figure  5- 
1.  The  spin  density,  per  carbon  atom,  on  rings  I (0.051)  and  III  (0.057)  has  been 
estimated  theoretically  to  be  approximately  four  times  larger  than  the  density  on  rings  II 
(0.026)  and  IV  (0.014)  [56].  This  concentration  of  spin  density  on  the  opposing  rings 
should  influence  the  geometry  of  the  g-tensor.  However,  the  influence  of  the  heavier 
nitrogen  atoms,  Nh  and  Niv,  may  outweigh  the  relatively  light  atoms,  with  small  spin- 
orbit  coupling  constants,  such  as  hydrogen  and  carbon.  The  spatial  spin  density  plots 
shown  in  figure  5-2  were  generated  by  Dr.  Alex  Angerhofer  from  DFT  calculations  to 
help  visualize  the  atoms  with  a higher  concentration  of  spin  density.  Figure  5-2(a) 
contains  a spatial  distribution  plot  of  the  alpha  spin  HOMO  state  (neutral  LUMO),  while 
the  total  spin  density,  which  is  the  difference  between  the  alpha  and  the  beta  spin 
densities,  is  shown  in  figure  5-2(b)  and  5-2(c).  In  figure  5-2(b)  and  5-2(c),  the  blue 
regions  of  highly  concentrated  spin  density  indicate  the  atoms  that  can  contribute  orbital 
angular  momentum  to  the  electronic  wavefunction.  Of  course,  atoms  such  as  nitrogen 
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and  oxygen  are  expected  to  have  a stronger  effect  on  the  g-tensor  than  carbon  and 
hydrogen.  The  nitrogen  atoms,  Nn  and  N|V,  are  indicated  with  the  arrows  in  figure  5-2(b) 
and  5-2(c).  According  to  the  DFT  predictions,  these  two  atoms  have  considerable  spin 
density,  which  should  lead  to  spin  orbit  interactions  in  this  dimension  of  g.  Lubitz  et  al. 
found  the  isotropic  hyperfme  couplings  to  be  AN(ii)  = 6.52  MHz  and  AN(iV)  = 5.86  MHz, 
with  negligible  couplings  on  N|  and  Nm,  using  TRIPLE  and  ENDOR  spectroscopy  in 
solution  [108].  However,  the  sheer  number  of  coupled  protons  and  the  weakly  coupled 
oxygen  atoms  on  ring  V and  the  carbon  2 acetyl,  could  introduce  more  spin-orbit 
interaction  along  the  Y dimension,  which  would  explain  the  predicted  g-tensor 
orientation  shown  in  figure  5-1 . A similar  g-tensor  orientation  has  been  determined 
through  experiment  and  theory  for  P700+  [124],  P865+  [33,73],  and  BChl  a [33]. 


*v 


Figure  5-2  BChl  a spin  density  plots  of  (a)  spatial  distribution  of  the  neutral  LUMO  (b) 
net  spin  density  spatial  plot,  (c)  average  net  spin  density  projection  with 
increasing  concentration  from  red  to  blue.  The  arrows  in  (b)  and  (c)  indicate 
nitrogens  II  and  IV,  which  should  strongly  contribute  to  the  g-factor. 

Isotropic  hyperfme  data  have  been  obtained  previously  with  EPR  [56],  ENDOR 
[56,108],  and  TRIPLE  [108]  resonance  techniques,  as  well  as  DFT  predictions  [109].  At 


X-Band  (9.5  GHz),  and  W-Band  (95  GHz),  frozen  solution  EPR  yields  a single 
featureless  line,  which  is  inhomogeneously  broadened  to  18  G (at  X-band)  by  unresolved 


hyperfme  couplings.  In  solution,  however,  this  anisotropic  interaction  is  averaged 
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beyond  the  timescale  of  the  experiment  and  the  isotropic  coupling  can  be  measured. 
Figure  4-1  contains  the  X-band  cw-EPR  spectrum  with  partially  resolved  hyperfine 
couplings.  These  splittings  represent  the  Fermi  contact  interaction  between  the  electron 
density  and  coupled  protons,  which  is  directly  proportional  to  the  electron  density  at  the 
carbon  nucleus.  The  spectrum  shown  in  figure  4-1  was  fitted  using  DOUBLET  software 
to  approximate  the  hyperfine  couplings.  The  values  obtained  were  in  agreement  with 
Fajer  et  al.,  who  performed  solution  state  EPR  and  ENDOR  experiments  on  selectively 
deuterated  BChl  a and  BPheo  a radical  anions  [56],  Fajer’s  work  concluded  that  the 
methyl  groups  of  carbons  1 and  5 produce  the  largest  couplings  of  3.1  G (8.7  MHz)  and 
3.3  G (9.2  MHz)  respectively  (see  figure  5-3  for  atomic  numbering).  These  couplings  are 
observed  in  solution  state  EPR  as  the  average  of  the  splittings  from  the  protons  on  each 
methyl  group.  The  methine  protons  (a,  P,  and  5)  produced  partially  resolved  ENDOR 
lines  with  hyperfine  couplings  of  2.0  G (5.6  MHz)  and  2.5  G (7.0  MHz)  with  relative 
intensities  1:2,  indicating  that  two  of  the  methine  protons  show  similar  hyperfine 
interaction  at  5.6  MHz,  and  one  shows  a stronger  interaction  at  7.0  MHz.  ENDOR  lines 
were  also  observed  at  0.6  G (1.7  MHz),  assigned  to  P-protons  (positions  3,  4,  7,  8),  and 
0.1  G (0.3  MHz),  which  were  assigned  to  the  remaining  weakly  coupled  protons.  The  X- 
band  EPR  experiment  shown  in  figure  4-1  identifies  the  radical  as  BChl  a radical  anion, 
and  measures  the  isotropic  hyperfine  couplings  for  reference.  The  simulation  in  figure  4- 
1 employs  the  following  couplings  to  fit  the  experimental  data:  Amethine=  3.0  G (8.4 
MHz),  Amethyi  = 2.4  G (6.72  MHz),  and  AN  = 2.2  G (6.16  MHz).  These  numbers  are  in 
agreement  with  published  data  [40],  and  are  presented  in  table  4-1. 
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Figure  5-3  Molecular  structure  and  atomic  numbering  of  BChl  a.  The  phytyl  tail 
(originating  from  carbon  7)  has  been  truncated  to  a methyl  group. 

In  order  to  build  on  the  isotropic  hyperfine  coupling  estimations  from  X-band, 
high  field  ENDOR  experiments  were  done  to  estimate  the  total  hyperfine  interaction. 
This  work  shows  orientationally  selected  ENDOR,  which  probes  the  total  (anisotropic 
and  isotropic)  hyperfine  coupling.  The  ENDOR  spectra,  shown  in  figure  4-8,  imply  that 
the  hyperfine  tensors  are  oriented  perpendicular  to  the  g-tensor  reference  frame 
[86,88,1 10],  This  relative  orientation  is  assumed  since  the  Az  component  of  the  A-tensor 
is  expected  to  lie  along  the  a bond  that  links  the  proton  to  the  conjugated  system.  As 
seen  in  figure  5-1,  the  g-tensor  maintains  gz  normal  to  the  plane  of  conjugation,  while  the 
protons  lie  within  the  gxy.  Also,  the  symmetric  spacing  of  the  ENDOR  transitions  about 
the  nuclear  Larmor  frequency  indicates  that  the  hyperfine  tensors  are  oriented  with  axes 
coincident  with  the  g-frame  [86,87],  In  this  radical,  hydrogen  and  nitrogen  nuclei 
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directly  interact  with  the  electronic  wavefunction,  but  only  proton  couplings  were 
measured  here,  due  to  experimental  limitations  on  the  spectrometer.  The  ENDOR  coil  of 
the  probe,  shown  in  figure  3-6,  was  optimized  for  protons.  Hydrogen  nuclei  are  about  10 

Y 

times  more  sensitive  to  NMR  since,  — = 10 . Due  to  the  difference  in  gyromagnetic 

Yn 

ratio,  nitrogen  nuclei  require  a much  larger  B2  field  to  induce  a nuclear  transition.  Since 
BChl  a radical  anion  has  strongly  and  weakly  coupled  protons,  a range  of  hyperfine 
interactions  is  observed. 

The  distant  or  weakly  coupled  nuclei  show  somewhat  isotropic  splittings.  That  is, 
the  hyperfine  interaction  is  weaker  in  general,  due  to  smaller  electron  density  at  the 
nuclei  (which  reduces  the  AjS0)  and  a larger  distance  from  the  electronic  wavefunction 
(which  reduces  Adjp).  The  anisotropic  (dipolar)  hyperfine  interaction  dies  off  as  r 3, 
where  r represents  the  averaged  distance  between  the  electron  density  and  the  coupled 
nucleus.  This  separation  dependence  is  shown  in  equations  2.10  and  2.11.  As  seen  in 
figure  4-8,  these  distant  protons  show  the  smallest  splittings,  between  0.20  MHz  and  1.0 
MHz.  It  is  reasonable  to  assume  that  the  ENDOR  intensity  at  the  Aj  = 0.20  MHz 
splitting  could  result  from  many  distant  protons  whose  coupling  constants  are  too  similar 
to  resolve  (see  figure  5-4).  These  transitions  could  result  from  protons  at  positions  3,  4, 

7,  8,  and  10,  which  have  been  estimated  experimentally  and  theoretically  to  have 
isotropic  couplings  between  0.50-1.8  MHz  [56,109].  The  presented  theoretical  DFT 
predictions  estimate  isotropic  couplings  of  1-2  MHz  for  these  protons.  The  splitting  of 
the  second  set  of  ENDOR  lines  increases  as  the  field  is  increased  from  gx  to  gz.  At  gx  the 
splitting  is  around  0.4  MHz  that  seems  to  increase  to  0.85  MHz  at  gz.  This  coupling  was 
fit  to  the  following  hyperfine  tensor:  A2  = [0.4,  0.85,  0.4]  MHz,  which  gives  A2(iSO)  = 0.55 
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MHz.  This  coupling  could  result  from  the  carbon  2 methyl  group,  which  is  known  to 
have  A2methyl(1S0)  = 0.4-0.5  MHz  [108].  The  third  coupling  is  weak,  and  becomes 
broadened  between  the  low  and  high  field  ENDOR  spectra.  The  approximate  splitting  of 
this  transition  ranges  from  3.0  MHz  at  gx  to  4.0  MHz  at  gz  and  was  fit  to  an  anisotropic 
hyperfine  tensor,  A3  = [3.0,  3.0,  4.0]  (in  MHz).  The  broadening  of  these  lines  through  the 
other  ENDOR  spectra  could  indicate  the  presence  of  another  pair  of  ENDOR  lines 
crossing  the  A3  lines.  So,  there  is  some  uncertainty  in  the  assignment  of  A3  to  this 
feature,  since  the  lines  cannot  be  distinguished  in  many  of  the  spectra.  The  fourth  pair  of 
ENDOR  transitions  has  a splitting  of  ~6  MHz  at  gx,  which  increases  to  ~7  MHz  at  gz. 
These  ENDOR  lines  fall  on  the  outer  wings  of  the  spectra,  which  appear  to  be  a 
superposition  of  two  or  three  unresolved  transitions.  The  expected  hyperfine  couplings 
for  the  methyl  and  methine  protons  are  similar  enough  to  fall  in  these  outer  wings.  So,  it 
is  conceivable  that  protons  from  a,  (3,  methyl  1,  and  methyl  5 are  superimposed  in  the 
region  of  A4  and  A5.  The  ENDOR  spectra  were  fitted  to  the  following  hyperfine  tensors: 
A4  = [6.6,  7.2,  9.2]  MHz,  where  A4(iS0)  = 7.7  MHz  and  A5  = [8.0,  9.2,  9.5]  MHz,  where 
Asuso)  = 8.9  MHz.  These  isotropic  couplings  have  been  compared  to  the  DFT  theoretical 
predictions  in  figure  5-5.  Since  ENDOR  experiments  are  blind  to  the  sign  of  a hyperfine 
coupling  constant,  these  numbers  could  be  positive  or  negative,  and  both  have  been 
plotted  in  figure  5-5.  The  best  agreement  with  the  literature  X-band  data  and  the 
theoretical  predictions  assign  A4  to  the  (3  proton  or  the  methyl  1 protons  (shown  in  table 
5-1  and  5-2).  The  isotropic  value  of  A5  agrees  with  the  a and  methyl  5 couplings. 
Measuring  the  ENDOR  at  30  K stopped  the  rotational  motion  of  the  methyl  groups, 
which  could  complicate  the  resolution  in  the  individual  protons’  resonances.  Higher 
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temperature  experiments  were  attempted,  but  the  EPR  intensity  was  too  weak  to  detect 
ENDOR  above  60  K.  Since  these  couplings  are  unresolved,  the  experimental  spectra  can 
only  be  estimated  through  speculation  and  theoretical  models. 


vrt  - v„  (MHz) 


Figure  5-4  BChl  a radical  anion  ENDOR  spectra  from  the  principal  g-tensor  components 
with  the  individual  transitions  labeled  A1-A5. 

Theoretical  (DFT)  methods  were  employed  to  investigate  the  hyperfine  couplings 
measured  in  the  hf-ENDOR  experiments.  Figure  5-5  plots  the  isotropic  hyperfine 
coupling  constants  calculated  with  the  indicated  basis  sets,  the  data  measured 
experimentally  by  Lubitz  et  al.  at  X-band  [108,109],  and  the  average  values  of  the 
hyperfine  tensors  derived  from  the  hf-ENDOR  experiments.  Agreement  between  the 
theoretical  predictions  and  the  X-band  experiments  is  good,  except  on  the  a and 
P protons.  The  theory  does  quite  well  for  the  methyl  protons,  but  deviates  from  the 
methine  protons:  cc  and  |3.  However,  the  DFT  calculations  presented  here  predict  the 
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correct  sign  of  the  hyperfine  coupling  constant,  which  is  an  improvement  over  pure  DFT 
functionals  [109]  as  well  as  semi-empirical  methods.  The  double  zeta  (DZ)  ADF 
calculations  show  the  best  agreement  for  protons  a,  p,  and  methyl  5,  while  the  best 
agreement  for  5 and  methyl  1 is  found  in  the  triple  zeta  (TZ)  calculation.  It  seems  that 
the  EPR  specific  basis  sets  and  higher  levels  of  theory  do  not  improve  the  predictions  in 
this  case.  The  predicted  anisotropic  hyperfine  tensors  are  given  in  tables  5-1  and  5-2, 
with  the  measured  experimental  data. 


Figure  5-5  Theoretical  and  experimental  isotropic  hyperfine  couplings  for  BChl  a radical 
anion.  The  hf-ENDOR  data  points  result  from  the  averaged  hyperfine  tensor 
components  of  A4  and  A5.  Since  the  sign  of  the  coupling  constant  is  uncertain 
in  ENDOR,  these  values  have  been  added  to  this  plot  with  both  sign 
possibilities. 
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Table  5-1  Isotropic  hyperfine  coupling  constants  estimated  with  DFT  and  X-Band  EPR 
for  BChl  a anion.  DFT  calculations  employed  a molecular  geometry  optimized 
with  one  THF  molecule  coordinated  to  Mg.  The  individual  component  of  each 
methyl  proton  tensor  has  been  averaged  to  anticipate  the  isotropic  nature  of  a 


Proton 

ADF 

ADF 

ADF 

Neese 

Neese 

Barone 

Exp. 

hf- 

(DZ) 

(TZ) 

(TZ, 

(DZ, 

(DZ, 

(EPRII) 

[109] 

ENDOR 

revPB) 

Ahlrich) 

IGLOII) 

a 

A iso 

-5.41 

-3.65 

-3.32 

-3.81 

-3.71 

-3.84 

-9.65 

-8.9? 

P 

A- 

*mso 

-5.43 

-3.65 

-3.08 

-4.75 

-4.64 

-4.76 

-6.91 

-7.7  ? 

5 

Ajso 

-9.55 

-6.81 

-5.52 

-10.72 

-10.25 

-10.45 

-6.23 

Methyl 

1 

A,av 

6.53 

6.64 

5.99 

6.26 

6.93 

7.39 

A2av 

7.01 

7.10 

6.06 

6.65 

7.34 

7.84 

A3av 

8.48 

9.35 

0.75 

9.24 

9.96 

10.36 

Ais„av 

7.60 

7.70 

6.49 

7.39 

8.08 

8.53 

7.63 

7.7  ? 

Methyl 

A,av 

7.84 

7.25 

6.86 

5.52 

6.16 

6.50 

A2av 

8.31 

7.66 

6.99 

5.91 

6.55 

6.93 

A3av 

9.82 

8.99 

8.33 

7.73 

8.39 

8.67 

A-  av 

^lSO 

8.81 

8.12 

7.14 

6.39 

7.03 

7.36 

9.19 

8.9  ? 
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Table  5-2  Hyperfine  tensor  data,  where  Ax  = Adjp  + AjS0 


Proton 

Position 

ADF 

(DZ) 

ADF 

(TZ) 

ADF 

(TZ, 

revPB) 

ADF 

(TZ, 

revPB 

-2p) 

Neese 

(DZ, 

Ahlrich) 

Neese 

(DZ, 

IGLOII) 

Barone 

(EPRII) 

hf- 

ENDOR 

a 

A-i 

-7.26 

-5.31 

-5.04 

-5.14 

-5.31 

-5.23 

-5.37 

A5 

8.0? 

^2 

-6.73 

-4.86 

-4.53 

-4.43 

-5.03 

-4.95 

-5.09 

9.2? 

a3 

-2.25 

-0.78 

-0.39 

-0.50 

-1.10 

-0.96 

-1.06 

9.5  ? 

P 

Ai 

-7.37 

-5.32 

-4.73 

-4.85 

-6.63 

-6.55 

-6.67 

A4 

6.6? 

A2 

-6.51 

-1.50 

-3.93 

-3.83 

-5.59 

-5.48 

-5.61 

7.2? 

A3 

-2.42 

-1.13 

-0.57 

-0.66 

-2.02 

-1.90 

-2.01 

9.2  ? 

5 

A, 

-13.24 

-10.53 

-9.07 

-9.27 

-15.69 

-15.20 

-15.41 

a2 

-11.12 

-7.83 

-6.56 

-6.36 

-11.80 

-11.33 

-11.59 

a3 

-4.30 

-2.06 

-0.94 

-1.12 

-4.67 

-4.22 

-4.36 

Methyl 

A* 

1 

A-i 

11.88 

11.99 

10.29 

10.20 

11.47 

12.63 

13.44 

6.6? 

A2 

12.42 

12.41 

10.66 

10.63 

11.78 

13.00 

13.84 

7.2? 

a3 

14.69 

14.72 

12.91 

12.81 

14.40 

15.67 

16.35 

9.2  ? 

Ai 

7.94 

8.28 

6.89 

6.78 

7.79 

8.41 

9.03 

A2 

8.34 

8.68 

7.30 

7.27 

8.11 

8.75 

9.46 

a3 

10.52 

10.84 

9.39 

9.33 

10.60 

11.27 

11.84 

Ai 

-0.25 

-0.36 

0.79 

-0.80 

-0.46 

-0.25 

-0.31 

a2 

0.27 

0.22 

0.21 

-0.20 

0.05 

0.26 

0.23 

a3 

0.26 

2.50 

0.20 

1.95 

2.74 

2.95 

2.89 

Methyl 

A5 

5 

Ai 

12.83 

11.49 

10.00 

9.88 

8.83 

9.74 

10.28 

8.0? 

A2 

13.34 

11.88 

10.40 

10.35 

9.23 

10.16 

10.75 

9.2  ? 

a3 

15.43 

13.68 

12.15 

12.09 

11.00 

11.96 

12.41 

9.5  ? 

Aj 

11.28 

10.93 

9.55 

9.42 

8.39 

9.27 

9.79 

A2 

11.72 

11.35 

9.98 

9.93 

8.82 

9.71 

10.29 

A3 

13.86 

13.12 

11.70 

11.64 

10.56 

11.48 

11.92 

Ai 

-0.60 

-0.67 

1.03 

-1.02 

-0.65 

-0.54 

-0.59 

a2 

-0.12 

-0.26 

0.59 

-0.6 

-0.33 

-0.22 

-0.27 

A3 

0.16 

0.16 

1.13 

1.13 

1.64 

1.73 

1.67 
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Bacteriopheophytin  a Radical  Anion 

The  high  field  EPR  spectrum  for  the  BPheo  a radical  anion  is  shown  in  figure  4- 
15,  and  the  g-tensor  data  of  repeated  experiments  are  tabulated  in  table  4-4.  The  sample 
for  each  experiment  was  prepared  in  the  same  way,  as  described  in  chapter  3,  but  may 
vary  in  concentration  since  the  yield  of  the  reduction  reaction  is  unknown.  The  g-values 
for  this  radical  provide  a unique  distinction  between  the  BPheo  a and  BChl  a radical 
anions  since  the  former  has  a more  anisotropic  g-tensor  than  the  latter.  BPheo  a radical 
anion  has  gz  = 2.002299  + 5xl0'5,  which  is  within  experimental  error  of  the  free  electron 
g-value,  ge  = 2.002319.  The  proximity  of  gz  to  ge  leads  to  the  assumption  that  gz  is  close 
to  perpendicular  to  the  conjugation  plane  (as  mentioned  in  the  BChl  a radical  anion 
discussion).  Figure  5-6  shows  the  proposed  orientation  of  g on  the  BPheo  a structure. 


Figure  5-6  Molecular  structure  and  atomic  numbering  of  BPheo  a with  the  estimated  g- 
tensor  orientation  for  BPheo  a radical  anion. 
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In  order  to  fit  the  hf-EPR  spectra,  an  anisotropic  linewidth  was  used,  which 
required  a broader  linewidth  for  gz  than  for  gx  and  gy.  This  result  could  originate  from 
anisotropic  hyperfine  interactions,  or  from  g-strain.  Since  the  gz  orientation  is  more 
likely  to  have  a distribution  of  solvent  and  intermolecular  interactions,  the  gz  feature  may 
become  broadened,  relative  to  the  other  components,  as  a result  of  a distribution  of 
environments  in  the  Z coordinate.  Such  interactions  have  been  observed  in  bacterial 
reaction  centers  and  quinone  anion  radicals  in  frozen  solution  and  in  single  crystals  [111- 
113].  The  possible  environmental  affects  on  the  gz  component  will  be  influenced  by  the 
nature  of  the  excited  states  of  the  radical.  That  is,  since  the  unpaired  electron  density  is 
located  in  the  pz  orbitals  of  the  carbon  atoms,  it  is  conceivable  that  environmental 
interactions,  with  other  BPheo  a radicals  or  solvent  (THF)  molecules,  could  broaden  gz. 
These  interactions  may  be  responsible  for  the  poor  fits  along  gz  on  the  hf-EPR  spectra. 
This  larger  uncertainty  in  the  simulations  translates  directly  to  the  uncertainty  of  the  g- 
tensor. 

The  values  of  gx  and  gy  were  measured  to  be  2.004459  + 5x10 5 and  2.003249  + 
5xl0'5  respectively.  Spin-orbit  interactions  lead  to  the  deviation  of  gx  and  gy  from  gz,  and 
this  interaction  is  heavily  weighted  by  the  spin-orbit  coupling  constant  (X)  of  the  nitrogen 
nuclei  of  rings  II  and  IV.  The  net  spin  density  on  Nn  and  Nrv  are  indicated  in  figure  5-7 
with  arrows.  Figure  5-7(a)  is  a plot  of  the  spatial  distribution  of  the  neutral  LUMO 
function  of  BPheo  a.  The  molecular  orientation  in  this  figure  is  identical  to  figure  5-6. 
Figures  5-7(b)  and  5-7(c)  were  generated  by  subtracting  the  alpha  spin  density  from  the 
beta  spin  density,  which  results  in  a net  spin  density  plot,  revealing  only  the  nuclei  with  a 
high  concentration  of  unpaired  spin  density.  This  net  spin  density  determines  the  g- 
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tensor  orientation.  From  figure  5-2  the  spin  density  appears  to  be  concentrated  along  the 

Y axis  of  the  molecular  frame.  This  prediction  resembles  the  BChl  a anion  predictions, 
which  is  believed  to  have  a similar  orientation  of  the  g-tensor  with  gx  along  the  molecular 

Y axis  between  rings  I and  III. 


(a) 

<< 

(b)  net  spin  density  (difference  between  the  alpha  and  beta  HOMO),  (c) 
projection  of  the  average  net  spin  density  distribution  with  increasing 
concentration  from  blue  to  red. 


The  BPheo  a radical  anion  shows  similar  orientationally  selected  ENDOR  spectra 
to  BChl  a,  with  subtle  differences.  This  radical  has  been  studied  a bit  more  closely  with 
theory  and  anisotropic  hyperfine  couplings  have  been  predicted  by  O’Malley  with  the 


B3LYP  level  of  theory  [47,1 14].  The  atomic  numbering  for  this  radical  is  shown  on  the 
structure  in  figure  5-8.  The  presented  data  show  superscripts  to  indicate  the  literature 
values,  and  subscripts  indexing  the  ENDOR  transition  number  from  the  experimental 
spectra.  The  smallest  splitting  was  fitted  to  a hyperfine  tensor,  A,  = [0.5,  0.5,  0.7]  MHz 
(see  figure  5-8  for  the  labeled  hyperfine  couplings).  This  tensor  leads  to  an  isotropic 
hyperfine  coupling  value  of  Ai(jS0)  = 0.57  MHz.  O’Malley  predicts  the  proton  of  the  ring 
I nitrogen  to  give  such  a coupling,  ANH(iso)  = 0.50  MHz,  with  a total  hyperfine  tensor,  ANH 
= [2.8,  0.2,  -1.5]  MHz.  The  prediction  of  such  a tensor  does  not  necessarily  conflict  with 
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Figure  5-8  Molecular  structure  and  atomic  numbering  of  BPheo  a. 
the  presented  ENDOR  data,  but  the  superposition  of  many  resonances  with  effective 
splittings  between  0.5  MHz  and  2.2  MHz  complicate  tensor  assignments.  The  second  set 
of  ENDOR  transitions  are  described  by  A2  = [0.7,  1.8,  1.8]  MHz  and  A2(iS0)  = 1.43  MHz. 
The  features  of  this  coupling  likely  result  from  the  proton  on  carbon  8,  which  O’Malley 
predicts  to  have  A8  = [-0.2,  -1.3,  -1.8]  MHz,  with  A8(lso)  = -1.1  MHz.  Again,  the 
uncertainty  in  the  experimental  data  increases  at  small  hyperfine  couplings,  which  may 
explain  the  difference  in  the  first  component  of  A8  and  A2.  The  orientation  of  this  tensor, 
relative  to  the  g-tensor  axes,  may  be  assumed  from  the  experimental  ENDOR  data.  The 
A splitting  appears  to  be  at  a maximum  at  the  edges  of  the  powder  pattern,  at  gx  and  gz, 
and  at  a minimum  at  gy.  This  would  indicate,  assuming  an  axial  tensor  with  Ax  < Ay  = Az, 
that  the  Ax  axis  lies  along  gy.  This  assumption  puts  Ay  and  Az  in  the  gxz  plane. 

The  third  set  of  resolved  ENDOR  lines  show  somewhat  isotropic  behavior,  with 
A3  = [2.0,  2.2,  2.0]  MHz,  which  gives  A3(iS0)  = 2.1  MHz.  O’Malley  predicts  similar 
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behavior  for  the  proton  on  carbon  3,  giving  A3  = [-1.3,  -2.2,  -2.5]  and  Al(lS0)  = -2.0  MHz. 
This  proton  is  bound  to  the  ring  II  pyrrole  and  is  not  predicted  to  have  much  spin  density 
[47,1 13].  However,  the  possibility  exists  that  this  feature  in  the  ENDOR  spectra  is  a 
superposition  of  several  resonances  from  similar  protons,  such  as  those  at  positions  3,  4, 
and  7.  These  protons  are  positioned  on  opposing  sides  of  the  molecule  (rings  II  and  IV), 
which  could  mean  that  their  A-tensors  are  approximately  180°  apart  in  their  orientation 
relative  to  the  g-frame.  This  assumption  is  made  since  the  z-component  of  the  A-tensor 
usually  lies  on  the  bond  between  the  proton  and  the  molecule,  where  the  largest 
component  of  the  A-tensor  should  result.  O’Malley’s  predictions  for  the  position  3,  4, 
and  7 protons  show  similar  anisotropies  and  similar  isotropic  couplings  with  an  average 
of  AiSo(ave)  = 2.5  MHz  + 0.5  MHz.  The  fourth  and  fifth  couplings  comprise  the  wings  of 
the  ENDOR  spectra  and  show  anisotropic  behavior,  with  A4  = [8.5,  7.8,  7.5]  MHz,  and 
A5  = [9.1,  8.8,  8.3],  with  A4(iso)  = 7.93  MHz  and  A5(iSO)  = 8.73  MHz.  These  resonances 
are  most  likely  due  to  the  methyl  protons  on  positions  1 and  5.  The  X-band  data  predict 
an  average,  isotropic,  methyl  contribution  of  8.4  MHz.  Thus,  it  is  plausible  that  the  outer 
wings  result  from  methyl  proton  resonance,  since  the  high-field  anisotropic  values 
average  to  8.33  MHz.  This  agreement  is  also  in  line  with  previous  work  at  X-band, 
which  similarly  finds  a value  of  8.4  MHz  [56],  Theoretical  predictions  for  the  methyl 
tensors  seem  to  be  much  more  anisotropic.  O’Malley’s  predictions  assume  an  average 
value  for  the  methyl  protons,  which  should  not  be  the  case  here,  since  the  ENDOR  data 
were  collected  at  10  K.  It  is  possible  that  the  resolution  of  the  couplings  for  these  protons 
could  be  improved  by  measuring  at  higher  temperatures.  Above  100  K the  methyl 
rotation  should  sufficiently  average  out  the  individual  proton  couplings. 
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Beyond  the  outer  wings  in  the  ENDOR  spectra,  a 13-15  MHz  splitting  seems  to 
be  present  as  shoulders  on  the  wings,  around  A5,  and  was  fitted  to  A6  = [15,  15,  13]  MHz. 
This  region  was  not  examined  in  all  the  spectra,  but  the  feature  is  present  in  each  of  the 
wide  scans  shown  in  figure  4-20.  A hyperfine  coupling  of  this  magnitude  is  predicted  in 
O’Malley’s  calculations  for  the  methine  proton  P where  A^  = [-3.4,  -8.8,  -12.1]  MHz  and 
A^iso)  = -8.1  MHz  [47].  For  such  a hyperfine  tensor  to  exist,  the  ENDOR  transitions 
must  be  very  weak,  or  very  broad  since  they  are  not  obvious  in  the  hf-ENDOR 
experiments  presented  here.  A strongly  anisotropic  tensor  such  as  this  would  result  in 
ENDOR  intensity  between  the  outer  wings  (A4)  and  the  central  feature  of  the  ENDOR 
spectra  (A3)  in  figure  5-8.  The  lack  of  intensity  in  this  region,  in  every  ENDOR 
spectrum,  seems  to  indicate  that  such  an  anisotropy  does  not  exist,  or  is  too  broad  and 
weak  to  be  detected.  However,  the  A6  splitting  could  be  due  to  methyl  protons. 

Measurement  of  the  ENDOR  at  10  K freezes  out  the  methyl  group  rotation.  In 
this  case,  the  individual  protons  on  each  methyl  group  will  have  different  orientations 
(i.e.  different  orientations  of  their  A-tensor),  as  well  as  different  hyperfine  couplings. 

This  scenario  should  add  a pair  of  lines  to  the  ENDOR  spectra  for  each  different  proton, 
unless  the  lines  are  not  fully  resolved.  According  to  theoretical  estimations  for  BChl  a 
radical  anion,  the  closer  a methyl  proton  is  to  being  in  the  plane  with  the  conjugation 
system,  the  smaller  the  total  hyperfine  interaction  (shown  in  table  5-2).  This  behavior 
results  from  the  nature  of  the  spin  density  on  this  radical,  which  is  present  above  and 
below  the  plane  of  conjugated  carbon  atoms.  Since  the  isotropic  couplings  are  an 
average  of  all  the  protons  of  a given  methyl  group,  a negligible  value  for  one  of  these 
protons  results  in  above  average  values  for  the  other  two.  In  other  words,  if  one  of  the 
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methyl  protons  at  10  K has  a negligible  hyperfine  coupling,  the  other  two  protons  must 
have  a coupling  much  larger  than  8-9  MHz,  which  is  the  averaged  methyl  coupling.  So, 
the  13-15  MHz  splitting  could  be  due  to  methyl  protons  that  are  out  of  the  plane  on 
methyl  1 or  5. 


Vrf  - vH  (MHz) 

Figure  5-8  Hf-ENDOR  spectra  of  BPheo  a radical  anion  with  labeled  hyperfine 
couplings. 
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Table  5-3  Hyperfine  tensor  data  for  BChl  a and  BPheo  a radical  anions.  The  tensors 


were  obtained  by  fitting  the  experimental  ENDOR  spectra  using  GENDOR 
[107].  All  simulations  used  Euler  angle  (j)  = 90°. 
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Bacteriochlorophyll  a Radical  Cation 

BChl  a is  a versatile  component  of  the  bacterial  reaction  center  of  Rb.  sphaeroides 
and  has  been  the  subject  of  much  low-field  EPR  and  ENDOR  research  [33-36,115].  This 
molecule  helps  collect  solar  energy  by  it’s  presence  in  the  light  harvesting  complexes 
(LH1  and  LH2),  and  is  involved  in  the  primary  charge  separation;  the  primary  donor  (D) 
is  a dimer  of  BChl  a molecules,  and  the  first  intermediate  acceptor  in  the  primary  charge 
separation  is  likely  to  be  a BChl  a monomer.  Essentially,  the  initial  charge  separation 
involves  the  interaction  of  the  cationic  (HOMO)  and  anionic  (LUMO)  states  of  BChl  a 
molecules.  The  cationic  species  is  a direct  probe  of  the  HOMO  wavefunction,  which  is 
the  first  species  created  in  the  primary  charge  separation,  D+.  Although  the  primary 
donor  is  a dimer  of  BChl  a molecules,  previous  work  has  shown  that  the  electronic 
structure  of  the  dimer  is  highly  asymmetric,  and  the  g-values  are  too  similar  to 
distinguish  between  D+  and  BChl  a [33].  High-field  studies  on  the  primary  donor,  and 
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monomeric  BChl  a cations  have  produced  isotropic  hyperfine  coupling  information,  as 
well  as  g-tensor  values  (Agx.z  = 0.0011)  and  orientations  from  single  crystal  experiments 
[33,35,74],  In  addition  to  the  small  g-anisotropy  of  this  radical,  the  hf-EPR 
measurements  were  complicated  by  the  presence  of  g-strain.  This  is  a broadening 
phenomenon  that  is  proportional  to  the  external  field.  As  the  resonance  conditions  of  the 
experiment  were  increased  from  0.35  T/9.5  GHz  to  7 T/210  GHz,  g-strain  became  more 
obvious  and  eventually  broadened  the  powder  pattern  several  hundred  gauss  at  210  GHz. 
This  broadening  effect  could  be  circumvented  by  deuteration  of  the  radical,  but  since  this 
option  was  not  immediately  available,  lower  frequency  measurements  were  pursued. 

X-band  proton  and  nitrogen  hyperfine  couplings  have  been  published  by  others 
[116-119]  and  calculations  have  shown  good  agreement  with  experimental  values  [109], 
According  to  Sinnecker  et  al,  the  electron  density  is  concentrated  on  the  beta  protons  of 
rings  II  and  IV  at  positions  3,  4,  7,  and  8 [109],  This  finding  is  in  contrast  to  the  radical 
anion’s  spin  density  distribution,  which  is  heavily  localized  on  the  alpha  protons  and  the 
position  1 and  5 methyl  groups.  For  the  cation,  these  beta  protons  show  isotropic 
hyperfine  couplings  of  1 1-16  MHz,  while  the  alpha  protons  have  much  smaller  couplings 
of  1-2  MHz  [67,136].  Figure  5-9,  contains  plots  of  the  spin  density  predicted  with  DFT 
calculations.  The  alpha  HOMO  level  is  shown  in  figure  5-9(a),  while  the  more  revealing 
plots  for  the  hyperfine  and  g-tensor  behavior  are  shown  in  5-9(b)  and  5-9(c).  These 
images  depict  the  predicted  net  spin  density,  which  explains  the  hyperfine  interactions 
and  the  origin  of  the  g-tensor.  Relative  to  BChl  a radical  anion,  the  cation  shows  a much 
more  even  distribution  of  spin  density  around  all  the  pyrrole  rings,  especially  on  the 
nitrogen  atoms.  This  observation  could,  in  principle,  explain  the  smaller  g-anisotropy 
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Figure  5-9  Spatial  spin  density  plots  generated  with  DFT  predictions  for  BChl  a radical 
cation,  (a)  The  alpha  HOMO  wavefunction,  (b)  depicts  the  net  spin  density, 
and  (c)  is  the  average  net  spin  density  plot  showing  the  increasing  spin 
concentration  of  blue  to  red. 

observed  for  the  cation.  The  spin  density  appears  to  be  similarly  localized  on  each 
of  the  four  pyrrole  nitrogen  atoms.  This  prediction  has  been  experimentally  justified  with 
the  hyperfine  coupling  measurements  of  the  nitrogens,  AN  = 2.7  MHz  + 0.5  MHz  [108], 
This  even  distribution  of  spin  density  should  influence  the  deviation  of  gx  and  gy  from  gz 
similarly.  The  carbonyl  group  on  ring  V also  appears  to  have  a substantial  amount  of  net 
spin  density,  which  adds  to  the  spin-orbit  coupling  potential  of  the  X dimension.  But, 
since  Nn  and  Niv  are  directly  linked  to  the  conjugation  system,  the  Y dimension  may 
have  more  spin-orbit  interaction  due  to  nitrogen.  The  spin  density  for  this  radical  seems 
to  be  concentrated  between  rings  II  and  IV.  High  field  studies  have  assigned  a g-tensor 
orientation  which  aligns  gx  with  X and  gy  with  Y [33],  The  BChl  a radical  anion,  on  the 
other  hand,  has  spin  density  concentrated  on  rings  I and  III,  which  should  orient  the 
anion’s  g-tensor  -90°  to  the  cation’s  g-tensor. 


The  W-band  echo  detected  field  sweep  of  BChl  a radical  cation  in  frozen  solution 
shows  little  Zeeman  resolution,  indicating  higher  fields  are  necessary  to  obtain  full 
orientation  selection  (shown  in  figure  4-23).  However,  the  possibility  of  partial 
orientational  selection  was  not  abandoned  and  Davies  and  Mims  ENDOR  spectra  were 
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taken  through  the  EPR  line.  Within  experimental  error,  these  ENDOR  spectra  show  no 
anisotropic  behavior.  This  result  is  not  surprising,  given  the  small  g-anisotropy  for  this 
radical,  Agx_z  = 0.0011  [33],  The  beta  protons,  at  positions  3,  4,  7,  and  8 show  up  weakly 
on  the  Davies  spectra  as  wings  at  As  14  MHz.  The  highly  dipolar  nature  of  this 
hyperfine  interaction  leads  to  very  broad  resonances  which  are  unresolved.  These 
ENDOR  frequencies  are  also  likely  to  be  slightly  different,  which  contributes  to  the  lack 
of  resolution.  The  position  1 and  5 methyl  resonances  are  expected  to  be  around  4.93 
MHz  and  9.62  MHz,  while  the  (3  and  8 proton  contributions  are  expected  around  1.30 
MHz  [67,109].  These  numbers  agree  very  well  with  the  measured  splittings  of  1.15 
MHz,  5.0  MHz,  and  9.5  MHz.  The  invariance  of  these  hyperfine  couplings  between  the 
different  ENDOR  spectra  indicates  that  the  g-tensor  resolution  was  insufficient  to  select 
orientations  from  the  powder  average. 

So,  the  BChl  a cation  and  anion  seem  to  have  very  different  spin  density 
distributions,  which  lead  to  their  very  different  spectroscopic  behavior.  The  pulsed 
ENDOR  experiments  presented  here  confirm  the  known  hyperfine  data  and  contribute  the 
knowledge  that  W-band  provides  insufficient  Zeeman  resolution  for  orientational 
selection.  In  order  to  conduct  orientational  selection  experiments  above  W-band, 
deuterated  samples  may  be  needed. 

Bacteriopheophytin  a Radical  Cation 

X-band  experiments  performed  on  this  cation  yield  an  unresolved  EPR  line, 
which  may  be  partially  resolved  at  W-band,  see  figure  4-27.  The  W-band  EPR  line  is 
broader  than  the  X-band  line,  and  is  beginning  to  show  some  asymmetry  on  the  low  field 
side,  possibly  indicating  spectral  separation  in  the  Zeeman  orientations.  However,  the 
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sharp  feature  at  33668  G may  result  from  another  radical  species  formed  during  the 
oxidation.  The  inaccuracy  of  the  magnetic  field  at  W-band  prevents  g-value 
determination  of  the  features,  which  could  help  identify  the  radical.  The  similarity 
between  the  pulsed  ENDOR  spectra  at  each  of  these  features,  shown  in  figures  4-28  and 
4-29,  indicates  the  presence  of  only  one  type  of  radical.  This  cation  has  been  partially 
characterized  by  J.  Fajer  at  X-band  using  deuteration  and  theoretical  techniques  [120]. 
Fajer  et  al.  report  gjso  = 2.0025  for  BPheo  a+,  which  is  in  agreement  with  this  work,  which 
finds  giso  = 2.00253  + 4x1 0’5  at  X-band.  The  measured  hyperfine  couplings  from  W-band 
ENDOR  experiments  are  presented  in  table  4-3.  It  has  been  predicted  with  theory  that 
rings  II  and  IV  carry  most  of  the  electron  spin  density,  which  leads  to  the  largest 
hyperfine  interactions  [120].  However,  the  hyperfine  couplings  from  the  protons  on  these 
rings,  1 1-20  MHz,  exceed  the  couplings  measured  here  [120],  So  the  splittings  observed 
in  the  W-band  ENDOR  experiments  result  from  the  methine  and  methyl  protons. 

The  lack  of  orientational  selection  in  these  experiments  prevents  the 
determination  of  hyperfine  tensor  components,  but  the  isotropic  constants  can  be 
assigned.  In  light  of  the  known  assignments  for  BChl  a radical  cation,  The  1.6- 1.3  MHz 
coupling  is  assigned  to  the  methine  protons,  [3  and  8.  These  protons  produce  1.3  MHz 
coupling  in  BChl  a radical  cation  [108,109].  The  3. 1-3.6  MHz  coupling  could  be  due  to 
the  remaining  methine  proton,  a,  which  shows  a 2.35  MHz  coupling  in  BChl  a radical 
cation.  This  larger  coupling  is  plausible  considering  the  absence  of  the  metal  atom 
should  distribute  more  electron  density  on  the  conjugated  atoms.  The  5.2  MHz  coupling 
closely  resembles  the  4.93  MHz  coupling  for  BChl  a radical  cation,  which  has  been 
assigned  to  the  position  1 methyl  group.  The  position  5 methyl  group  shows  a coupling 
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of  9.62  MHz  for  BChl  a cation,  so  the  7. 7-8.7  MHz  splitting  is  assumed  to  belong  to  the 
protons  of  this  group. 

High  field  EPR  measurements  on  BPheo  a radical  cation  were  attempted  with 
similar  results  to  BChl  a radical  cation.  The  presence  of  g-strain  prevented  full  resolution 
in  the  g-tensor,  and  thus,  orientationally  selected  ENDOR  experiments.  The  spectra 
measured  here  show  that  W-band  does  not  provide  sufficient  Zeeman  resolution  to  select 
individual  orientations  from  the  powder  average.  Deuteration  could  possibly  improve  on 
the  information  contained  in  W-band  spectra,  as  was  suggested  for  BChl  a cation.  But 
the  measured  isotropic  couplings  indicate  a similar  spin  density  distribution  to  BChl  a 
radical  cation,  with  slightly  more  density  on  the  position  1 methyl  and  a slightly  smaller 
density  on  the  position  5 methyl.  Also,  the  methine  proton,  a,  seems  to  be  more  strongly 
coupled  in  BPheo  a cation  than  in  BChl  a cation. 

Chlorophyll  Radical  Anions 

Fully  resolved  g-tensor  measurements  on  the  radical  anions  of  Chi  a and  Chi  b 
have  been  taken  and  are  shown  in  figures  4-30  and  4-32,  and  the  g-tensor  data  is 
tabulated  in  table  4-4.  The  magnitude  and  anisotropy  of  the  g-tensors  are  similar  between 
these  two  chlorophyll  species.  The  g-tensor  of  doublet  radicals  is  known  to  be  very 
sensitive  to  the  molecular  structure  of  the  radical  [74,80,121,122],  and  intermolecular 
interactions  [123].  The  structural  difference  between  these  two  chlorophylls  amounts  to  a 
methyl  group  on  carbon  3 for  Chi  a,  where  Chi  b has  an  acetyl  group.  Zeeman  resolution 
was  achieved  by  330  GHz,  which  seems  to  indicate  the  lack  of  g-strain,  which  has  been 
demonstrated  by  the  Chi  a+  [74,124].  Since  the  high  field  EPR  experiments  were  done 
on  frozen  solutions,  the  principal  g-tensor  values  were  measured,  but  assignment  of  an 
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orientation  relative  to  the  molecular  frame  requires  additional  methods.  However,  it  is 
possible  to  make  certain  conclusions  from  the  g-tensor  data. 

The  proximity  of  the  chlorophyll  radical  anions’  gz  to  the  free  electron  value  is 
expected  since  these  species  are  conjugated  7t-radicals  [125].  The  Chi  a radical  anion 
showed  gz  = 2.002024  + 8xl0"5,  and  for  Chi  b radical  anion,  gz  = 2.001896  + 7xl0"5.  The 
values  of  gz  indicate  the  shift  of  the  entire  g-tensor  between  the  two  radicals.  This  shift  is 
expected  from  the  X-band  experiments,  where  AgiS0(Chl  a - Chi  b)  = 1.8xl0"4.  The 
anisotropy  between  the  principal  components  of  g is  similar  between  these  two  radicals, 
but  the  absolute  value  of  g is  shifted.  This  shift  in  the  high  field  g-values  (Agave(Chl  a - 
Chi  b ) = 1 . 14x  10"4)  exceeds  the  uncertainty  in  the  field  calibration  standard  ( + 5xlO'5 ) 
and  exceeds  the  expected  shift  of  g due  to  intermolecular  affects  with  THF  or  residual 
water  (on  the  order  of  lxlO’4  + 5xlO"5).  These  environmental  effects  would  most  likely 
affect  the  gz  component  more  than  the  other  components  of  g,  since  the  unpaired  electron 
density  occupies  this  dimension  generated  by  the  pz  atomic  orbitals. 

These  two  molecules  differ  in  structure  on  carbon  3,  where  Chi  a has  a methyl 
group  and  Chi  b has  an  acetyl  group.  If  these  functional  groups  have  considerable  spin 
density,  then  the  presence  and  orientation  of  the  acetyl  group  could  heavily  influence  the 
value  of  the  g-tensor.  The  presence  of  an  oxygen  atom  can  greatly  affect  the  value  and 
orientation  of  the  g-tensor  components,  as  has  been  shown  on  quinone  radicals  [ref]. 
However,  the  extent  of  the  affect  on  the  g-tensor  is  directly  related  to  the  amount  of  spin 
density  localized  on  the  oxygen,  which  appears  to  be  quite  small  for  the  chlorophyll 
radical  anions.  The  affect  of  the  acetyl  orientation  on  the  g-tensor  has  been  predicted 
using  the  BLYP  level  of  DFT  theory,  on  other  photosynthetic  ions,  and  can  shift  the  g- 


105 


tensor  components  by  up  to  17xl0‘4  [33,111].  The  results  presented  here  indicate  that 
this  functional  group  has  a small  affect  on  the  total  anisotropy  of  the  g-tensor,  since  the 
difference  in  the  total  anisotropies  is  given  by:  |Agx.z(Chla  - Chlb)|  = 9xl0'5.  This 
difference  in  the  g-anisotropy  is  on  the  order  of  the  error  in  the  g-tensor  measurements, 
which  is  5xl0'5  to  lxlO"4.  Thus,  the  interaction  of  the  acetyl  group  with  the  electron 
density  seems  to  be  negligible.  If  it  is  assumed  that  carbon  3 falls  on  the  Y-axis  of  the 
molecular  reference  frame,  the  g-tensor  component  closest  to  this  axis  should  be  the  most 
strongly  affected  by  the  presence  and  orientation  of  the  acetyl  group.  But  since  the  g- 
anisotropies  are  similar  between  the  Chi  a and  Chi  b radical  anions,  a g-tensor  orientation 
cannot  be  assumed  from  this  data. 

Density  functional  estimations  of  the  spin  density  on  the  Chi  a radical  anion 
predict  the  acetyl  group  to  have  negligible  spin  density  [126],  This  result  is  surprising 
considering  that  chlorin  pigments,  such  as  chlorophylls,  have  a conjugation  system  that 
extends  to  rings  I,  II,  and  III.  This  conjugation  brings  carbon  3 into  an  alpha  position  (on 
ring  II),  where  it  could  possess  considerable  electron  density.  However,  the  negligible 
affect  of  the  acetyl  group  on  the  g-tensor  agrees  with  the  experimental  data  and  the 
theoretical  predictions  of  Sinnecker  [126]  and  O’Malley  [127]  who  have  calculated  a 
negligible  spin  density  on  ring  II.  This  result  is  in  contrast  to  the  bacterial  anions,  which 
are  predicted  to  use  the  carbon  2 acetyl  group  in  the  primary  electron  transfer  events 
[47,114],  This  difference  may  indicate  an  evolution  of  a more  specialized  chlorophyll. 
That  is,  the  electron  density  and  chemically  active  groups  lie  along  the  X-axis 
(connecting  rings  I and  III).  In  this  case,  the  acetyl  group  on  carbon  3 may  not  be  directly 
involved  in  the  charge  separation  chemistry. 
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The  orientation  of  the  g-tensor  has  been  estimated  for  Chi  a cation  to  have  gx 
connecting  rings  I and  III,  gy  connecting  rings  II  and  IV,  and  gz  approximately  normal  to 
the  conjugation  plane  [33,74,1 13].  This  orientation  seems  to  follow  the  predicted  spin 
density  symmetry  on  the  Chi  a anion  species  as  well,  which  maintains  subtle  differences 
from  the  cationic  wavefunction,  but  maintains  a higher  concentration  of  spin  density  on 
rings  I and  III,  and  a relatively  lower  concentration  on  rings  II  and  IV  [121].  The 
orientation  of  the  g-tensor  has  been  shown  theoretically  to  be  very  sensitive  to  the 
geometry  of  the  radical  [41,80].  The  presence  of  g-strain  impaired  the  g-tensor 
measurement  of  the  cation,  which  led  to  the  study  of  deuterated  material  [74].  This  g- 
strain  is  a broadening  phenomenon  that  increases  the  linewidth  as  the  field  is  increased. 
The  g-anisotropy  for  the  Chi  a radical  anion  exceeds  that  of  its  cationic  counterpart, 
which  most  likely  indicates  a similar  trend  in  spin  density  distribution  to  the  bacterial 
ions.  This  distribution  entails  strongly  coupled  rings  I and  III,  with  a more  symmetric 
distribution  on  all  of  the  pyrrole  nitrogens,  Nj  - Njv. 

The  variation  in  experimental  g-tensor  data  for  Chi  a can  be  rationalized  through 
solvation  effects.  It  is  known  that  solvent  ligation  as  well  as  aggregation  between 
radicals  and  their  neutral  counterparts  can  affect  the  electronic  structure  of  the  radical 
species,  which  directly  affects  the  g-tensor  [79,1 12],  These  effects  are  rooted  in 
electronic  interactions  such  as  hydrogen  bonding  or  n-n  stacking  (which  are  purely 
dipolar  interactions).  This  additional  dipolar  interaction  can  account  for  shifts  in  the  g- 
value  on  the  order  of  10  4 [1 13].  The  presence  of  heavy  atoms  can  increase  the  shift  for 
such  an  interaction.  The  solvent  for  all  the  anion  studies  was  tetrahydrofuran  (C4H80), 
which  contains  an  oxygen  atom  that  can  serve  as  the  fifth  ligand  to  the  central  Mg  atom 
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in  the  chlorophylls.  Although  water  was  rigorously  removed  from  the  system  (as 
described  in  chapter  3),  the  presence  of  residual  water  in  the  sample  could  promote 
aggressive  hydrogen  bonding  and  varied  ligation  to  Mg.  This  interaction  was  observed  as 
early  as  the  1970’s  to  explain  discrepancies  in  the  isotropic  g-value  for  BChl  a radical 
anion  [36],  The  principal  g-tensor  components  have  a relatively  small  uncertainty  (7x10’ 

3 + 2xl0"5),  but  the  possibility  of  environmental  affects  resulting  from  the  nature  of  the 
sample  increase  this  value  (up  to  2xl0"4).  This  level  of  uncertainty  is  still  very  small 
relative  to  the  g-values  and  should  not  take  away  from  the  utility  of  these  data. 

The  primary  charge  separation  in  photosystem  I is  still  under  investigation  with 
optical  and  magnetic  resonance  experiments  to  study  the  radical  species  formed 
[128,129],  The  primary  acceptor,  A0,  has  been  determined  by  crystallography  to  be  a Chi 
a monomer  [130],  which  takes  on  an  electron  during  the  primary  charge  separation.  This 
charge  separated  state  is  very  short-lived,  making  it  difficult  to  study  in  vivo.  Knowledge 
of  the  g-values  and  anisotropies  of  the  Chi  a radical  anion  could  be  useful  in  one  of  two 
ways:  the  identity  of  an  in  vivo  radical  may  be  determined  by  its  g-tensor,  or  the  extent 
of  the  interaction  between  the  radical  and  the  environment  may  be  assessed.  This  g- 
tensor  data  on  the  individual  anion  species  should  allow  further  identification  of  the 
radical  formed  during  the  primary  charge  separation,  A0,  which  has  never  been 
successfully  measured  with  high  field  EPR.  This  data  also  provides  a reference  frame  for 
the  measurement  of  environmental  interactions  within  the  RC.  It  may  be  possible  to 
estimate  the  extent  of  protein  interactions  by  the  difference  between  the  in  vivo  and  in 
vitro  g-tensor  values.  High-field  EPR  [72]  and  DFT  [49,83]  methods  have  been  used  to 
investigate  protein  interactions  with  organic  radicals.  The  presence  of  hydrogens  on  the 
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chlorophyll  structure  allows  hydrogen-bonding  possibilities  to  the  surrounding  proteins, 
which  could  affect  the  energetics  of  the  radical  species. 

Due  to  the  relatively  small  g-tensor  anisotropy  (~10'4),  as  well  as  from  the 
predicted  atomic  spin  densities  [127],  that  the  oxygen  and  nitrogen  atoms  have  a 
relatively  small  concentration  of  electron  density.  The  uncertainty  in  the  gz  component 
seems  to  be  much  larger  than  the  uncertainty  in  gx  or  gy.  This  could  result  from  the 
dipolar  solvent  ligation  that  is  sterically  limited  to  the  direction  of  gz.  This  ligation 
should  be  prominent  in  the  Mg  containing  pigments,  and  thus  could  explain  the 
heightened  variation  of  the  gz  component  for  the  metallated  pigments.  The  pheophytin 
radicals  were  not  measured  here,  but  the  heightened  gz  broadening  is  present  in  the 
bacterial  pigments  between  BChl  a and  BPheo  a radical  anions.  It  must  be  disclosed  that 
ligation  along  gz  does  not  guarantee  that  this  component  will  be  most  effected  by  the 
interaction.  However,  since  the  spin  density  is  largely  localized  in  the  pz  orbitals,  gz  is 
likely  to  be  more  affected  by  electrostatic  interactions  along  this  dimension  than  gx  or  gy. 


CHAPTER  6 
CONCLUSIONS 


The  high  field  EPR  data  report  the  first  full  characterization  of  the  electronic 
Zeeman  interaction  for  the  radical  anions  of  BChl  a,  BPheo  a,  Chi  a,  and  Chi  b (shown  in 
table  6-1).  The  BChl  a and  BPheo  a radical  cations  showed  unresolved,  g-strain 
broadened  EPR  spectra  at  240  GHz,  which  limited  study  to  95  GHz.  At  this  frequency, 
no  orientation  selection  was  achieved,  which  prevented  g and  A-tensor  measurement. 

The  high  field  EPR  on  the  radical  anions  estimates  the  nature  of  the  electron  density  and 
also  firms  the  differences  between  the  radical  anion  species  studied.  The  bacterial 
pigments  show  that  the  presence  of  the  central  Mg  atom  affects  gave  (giSO)  more  than  the 
g-anisotropy,  Agx.z.  However,  Ag  is  a more  reliable  experimental  signature  for  any 
radical,  since  the  magnitude  of  giso  is  limited  by  the  quality  of  the  field  standard.  The 
orientation  of  the  g-tensor  appears  to  be  similar  between  all  the  radical  anions,  which 
have  an  electron  spin  density  localized  on  rings  II  and  IV.  This  result  is  not  surprising 
considering  the  species’  structural  similarities.  Also,  the  spin  density  on  the  carbon  and 
hydrogen  atoms,  rather  than  the  nitrogen  and  oxygen  atoms,  dictates  the  g-tensor 
orientation.  This  behavior  is  seen  in  the  radical  cations  as  well,  where  a much  more 
symmetric  distribution  of  spin  density  results  in  a much  smaller  g-anisotropy.  Further 
work  should  be  done  to  measure  the  radical  cations’  g-tensors,  and  ultimately  A-tensors 
with  high  field  ENDOR. 
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These  experiments  could  also  be  extended  to  study  Chi  a and  Chi  b,  and  possibly 
Pheo  a and  Pheo  b radical  anions  and  cations.  The  exclusive  presence  of  the  a 
derivatives  in  plant  reaction  centers  merits  the  investigation  of  the  differences  between 
the  a and  b derivatives.  In  light  of  the  g-anisotropy  observed  for  Chi  a and  Chi  b radical 
anions,  the  high  field  ENDOR  experiments  should  be  straightforward  and  similar  to  the 
procedures  presented  here  for  BChl  a and  BPheo  a radical  anions.  The  study  of  the 
cations  will  likely  be  less  straightforward  since  deuteration  would  enhance  the  g-tensor 
resolution,  but  would  complicate  ENDOR  measurements.  If  selective  deuteration  were 
available  and  precise,  this  technique  would  be  ideally  suited  for  these  experiments. 

Table  6-1  Radical  anion  g-values  measured  with  high  field  EPR.  These  values 


correspond  to  the  high  field  EPR  spectra  shown  in  chapter  4 and  are  the  most 
reliable  of  the  experimental  data  for  each  species.  Reliability  was  determined 
by  the  quality  of  the  fitting,  which  is  indicated  by  the  experimental  error. 


Sample 

8xx 

gyy 

gzz 

(gxx'gyy) 

(gyy'gzz) 

Microwave 

Frequency 

(GHz) 

Bacteriochlorophyll  a ' 

2.004257 3 

2.0031 50 a 

2.002513“ 

0.001107 

0.000637 

210.37 

Bacteriopheophytin  a ' 

2.004459 a 

2.003249“ 

2.002299“ 

0.001210 

0.000950 

326.13 

Chlorophyll  a ' 

2.004012“ 

2.0029 12 a 

2.002027 d 

0.001100 

0.000885 

326.11 

Chlorophyll  b ' 

2.003964 b 

2.002779 b 

2.001 896 c 

0.001185 

0.000883 

326.13 

The  superscripts  a-f  denote  the  experimental  error  in  that  measurement:  a = 5xl0'5,  b = 6xl0'5, 
c = 7xl0‘5,  d = 8xl0‘5. 


The  high  field  ENDOR  spectra  provide  the  first  experimental  measurement  of  the 
anisotropic  hyperfine  interactions  for  BChl  a and  BPheo  a radical  anions.  This  data 
provides  experimental  estimates  of  the  hyperfine  tensors  and  their  orientations  relative  to 
g.  The  resolution  of  these  ENDOR  spectra  is  limited  by  the  ENDOR  linewidth  and 
similarity  of  the  hyperfine  splittings.  Out  of  the  12-16  possible  couplings,  5-6  are 
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observable.  Granted,  7-9  of  these  expected  couplings  have  AjSO  < 2 MHz.  The  linewidths 
of  the  ENDOR  transitions  ranged  from  0.2  MHz  to  2 MHz.  Such  broadening  may  result 
from  the  distribution  of  molecular  orientations,  which  was  accounted  for  using  gaussian 
lineshapes  in  the  fitting  procedures.  Regardless,  the  high  field  ENDOR  spectra  present 
orientationally  selected  data,  which  provides  experimental  data  for  parameters  that  have 
only  been  theoretically  predicted  in  the  past.  Although  quantitation  of  all  hyperfine 
tensors  was  complicated  by  poor  resolution  between  transitions  occurring  at  similar 
ENDOR  frequencies,  estimations  of  the  tensor  components  and  trends  in  the  splitting 
were  obtained.  The  largest  splittings,  A4  and  A5  for  BChl  a,  and  A5  and  A6  for  BPheo  a , 
are  tentatively  assigned  to  the  methyl  protons  of  positions  1 and  5 (shown  in  table  6-2). 
These  tensor  predictions  may  be  inaccurate  due  to  the  overlap  of  several  lines.  For 
example,  the  outer  band  of  the  BChl  a ENDOR  spectra  (figure  5-4)  is  predicted  to  result 
from  two  couplings,  A4  and  A5.  The  BPheo  a radical  anion  is  expected  to  behave 
similarly,  but  this  outer  band  is  totally  unresolved,  thus  limiting  the  assignment  of  an  A- 
tensor  (figure  5-9).  The  methine  protons  show  similar  average  hyperfine  couplings,  but 
are  predicted  by  DFT  to  show  larger  anisotropy  than  the  methyl  protons.  This 
assignment  may  be  investigated  further  by  examining  the  temperature  dependence  of  the 
ENDOR  spectra.  This  rotation  would  average  out  the  individual  methyl  proton 
interactions,  which  should  sharpen  the  ENDOR  lines  resulting  from  methyl  protons. 

Also,  partial  deuteration  could  help  determine  the  nature  of  the  hyperfine  interactions  by 
removing  the  ENDOR  lines  of  the  isotopically  substituted  protons  from  the  spectra. 

The  presented  high  field  ENDOR  experiments  are  truly  groundbreaking  and 
demonstrate  the  possibility  of  measuring  hyperfine  tensors  in  frozen  solution,  despite  the 
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Table  6-2  Experimental  hyperfine  data  for  BChl  a and  BPheo  a radical  anions. 


Proton 

BChl  a 

Anion 

BPheo  a 

Anion 

X-band 

Data 

[109] 

DFT 

Calcs. 

[A.Angerhofer] 

hf-ENDOR 

Data 

X-band 

Data 

[109] 

DFT 

Calcs. 

[47] 

hf-ENDOR 

Data 

Methyl  1 

A, 



6.64 

6.6 



8.0 

8.5 

a2 

- 

7.31 

7.2 

- 

6.0 

7.8 

a3 

- 

9.35 

9.2 

- 

5.5 

7.5 

A- 

^lSO 

7.63 

7.70 

7.7 

7.10 

6.5 

7.9 

Methyl  5 

Ai 

— 

7.25 

8.0 

— 

12.0 

9.1 

A2 

- 

7.66 

9.2 

- 

9.7 

8.8 

a3 

- 

8.99 

9.5 

- 

9.1 

8.3 

Aiso 

9.19 

8.12 

8.9 

8.32 

10.2 

8.7 

small  Zeeman  and  hyperfine  anisotropy  of  doublet,  7t-radicals  with  relatively  weak  spin- 
orbit  and  hyperfine  interactions.  Knowledge  of  these  g and  A-tensors  should  be 
immediately  useful  for  in  vivo  studies.  The  radicals  formed  during  the  primary  charge 
separation  may  be  identified  by  their  electronic  signature,  which  has  been  approximated 
in  this  work  in  vitro.  The  protein  subunits  of  the  reaction  center  could  influence  the  in 
vitro  signature  through  electrostatic  interactions  such  as  H-bonding.  Such  interaction 
may  complicate  the  in  vivo  identification  of  the  radical  species,  but  at  the  same  time  may 
be  useful  in  studying  the  environmental  influence  on  the  electronic  structure.  Single- 
crystal measurements  may  provide  fully  resolved  hyperfine  tensor  information  since  the 
dominant  broadening  affect  in  these  ENDOR  spectra  appears  to  arise  from  the 
distribution  of  environments  and  geometries  of  the  radicals  in  frozen  solution.  Such 
affects  would  be  eliminated  in  a single-crystal  sample,  but  preparing  such  a sample  is 
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currently  not  routine.  However,  crystallization  of  the  entire  reaction  center  is  possible 
and  produces  a highly  ordered  system  due  to  the  protein  structure  [33,73], 


APPENDIX 

ABBREVIATIONS  AND  CONSTANTS 


Abbreviation 

Definition 

A 

Amn 

ADF 

Pe 

Pn 

Ba 

Bo 

B, 

b2 

BChl  a 
BPheo  a 
Chi  a 
Chi  b 
cw 
D 

DFT 

DOUBLET 

EMR 

ENDOR 

EPR 

ESR 

ESEEM 

EFT 

$A 

g 

gave 

ge 

giso 

gmn 

GENDOR 

h 

H 

hfc 

HOMO 

I 

L 

LUMO 

Hyperfine  coupling  tensor 

Component  of  the  hyperfine  coupling  tensor 

Amsterdam  Density  Functional  Theory 

Bohr  Magneton:  9.27402  x 10"24  JT'1 

Nuclear  Magneton:  5.0507866  x 10"27JT' 

Bacteriochlorophyll  molecule  on  A-branch  of  bacterial  RC 

Static  magnetic  field 

Magnetic  component  of  microwave  radiation 

Magnetic  component  of  radiofrequency  radiation 

Bacteriochlorophyll  a 

Bacteriopheophytin  a 

Chlorophyll  a 

Chlorophyll  b 

Continuous  Wave 

Primary  Donor 

Density  Functional  Theory 

Simulation  program  used  to  fit  the  experimental  EPR  data 

Electron  Magnetic  Resonance 

Electron  Nuclear  DOuble  Resonance 

Electron  Paramagnetic  Resonance 

Electron  Spin  Resonance 

Electron  Spin  Echo  Envelope  Modulation 

Fast  Fourier  Transform 

Bacteriopheophytin  molecule  on  A-branch  of  RC 
g- tensor 

Average  g-value:  (gx+gy+gz)/3 

Free  electron  g-value:  2.002319304386(20)  [76] 

Isotropic  g-value 
Component  of  the  g-tensor 

Simulation  program  used  to  fit  the  high  field  ENDOR  data 

Planck’s  constant:  6.6260755(40)  x 10"34Js 

Hamiltonian  operator 

Hyperfine  coupling  constant 

Highest  Occupied  Molecular  Orbital  (of  neutral) 

Nuclear  spin  moment 

Orbital  angular  momentum 

Lowest  Unoccupied  Molecular  Orbital  (of  neutral) 
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M 

Bulk  spin  magnetization 

Ms 

Projection  of  electron  spin  moment  on  quantization  axis 

MO 

Molecular  orbital 

vH 

Proton  Larmor  frequency 

Vrf 

RF  frequency 

NHMFL 

National  High  Magnetic  Field  Laboratory 

NMR 

Nuclear  Magnetic  Resonance 

Bb 

Bohr  magneton:  9.2740154(31)  x 10"24  JT"1 

PS  I 

Photosystem  I 

PS  II 

Photosystem  II 

Qa 

Quinone  acceptor  on  A-branch  of  bacterial  RC 

RC 

Photosynthetic  Reaction  Center 

S 

Electron  spin  angular  momentum 

UV-Vis 

Ultraviolet  - Visible  absorption  spectroscopy 

W-band 

Naming  convention  for  ~95  GHz 

X-band 

Naming  convention  for  ~9.5  GHz 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 
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